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ON PRIMITIVE PYTHAGOREAN TRIANGLES WITH EQUAL PERIMETERS 


Leon Bernstein* 


58, Sokolov Street, #5, Tel-Aviv, Israel 
(Submitted November 1983) 


Dedicated to my wtfe Pesta and my son John 


1. Introduction 


A triple (x, y, 2) of natural numbers is called a Pythagorean Triangle if 
x, y, 2 satisfy the Pythagorean equation 


x2 + y2 = 22, 


The triple (x, y, 2) is a Primitive Pythagorean Triangle (PPT) if x, y, 2 have 
no common factor greater than 1. If x is assumed to be odd, the set of PPT's 
can be generated by the set of pairs of natural numbers (u, v) satisfying 


WSR e105 Cis OY 21, 4-4 =. 1 Giod 2), €1) 
the well-known generating formula being 
(x, ys 2) = (u2 - v2, 2uv, u2 + v2). 


The pair (u, v) is called the generator of the PPT (x, y, 2). 


In terms of the generator, the perimeter S of (x, y, 2), S=x+yt+2, may 
be expressed as 


S = 2u(u + v). 


Denote by H the set of all such perimeters. Let H, be the subset of H defined 
by the relation: S € H, if S is the perimeter of exactly k PPT's. 

It is not difficult to show that H, is an infinite set, i.e., there is an 
infinite set of PPT's each one of which has a perimeter not shared by any other 
PPT. The surprising fact that Hj is also an infinite set is proved in [1]. It 
is the main purpose of this paper to prove that H, is an infinite set for any 
k, k = 33; see Proposition 3.3 below. The proof may appear to be const -uctive, 
but it is ultimately seen to depend on a known existential Theorem of imalytic 
number theory, the so-called modern version of Bertrand's postulate. 

Necessary conditions for the construction of k PPT's with equal pe cimeters 


are given in the next section. That the conditions can be met is shorn in the 
proof of Proposition 3.3. 


2. A Constructive Device 


Let us first construct k different generators (u, v) of PPT's with equal 
perimeters. 


*This paper is the final version of two papers submitted for publication by Leon Bernstein before 
he died on March 12, 1984, of a cerebral hemorrhage. It benefitted from the advice of a number of 
anonymous referees. 
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ON PRIMITIVE PYTHAGOREAN TRIANGLES WITH EQUAL PERIMETERS 


Proposition 2.1: Let B,, Bo,..., By be kK (k 23) odd positive integers, pairwise 
relatively prime, By < By <q sae < Bys and 


B, < B,v2. (2) 
Let 


k 
A, = [] 8; and u, = A,/B, for t€ T, T= {1, 2, ..05 k}. 


t=1 
Assume there exists an odd positive integer P, satisfying the two conditions 
Brag ip ree, (3) 


BoD aie By BYES vise hp) 
By ao By : “3 


If v, = P.B, - uz» t © T, then the pairs (u,, v,) are generators of k PPT's hav- 
ing equal perimeters S, S = 2P,A,. 


Proof: We show first that (u,, v;) is the generator of a PPT for each t ET, 


i.e., that (u,, v4) satisfies (1). From the definitions of u,, vz, it follows 


that 
Uj, >Up > sts >U, and Vy < Vy < tt" < Vy. (5) 


Since by (4), 
DV, = BB, - u, = BB, - BBs «++ By > 0; 


it follows from (5) that v; > 0 for t € T. Moreover, it follows from (5) that 
Uz, > Ves EET, provided u, > v,. And this is a consequence of (4): 
yy +o - 2u, - PB, > (2A,/B,) - 2B,B, «-- By} = 0. 
Thus, be > Vis CET. 
Next, (uz, Ve) = 1 if and only if (uz, uz + Vt) = (Ay/Be, Py Be) = 1, which 


is true since, by assumption, (u;, P;,) = 1 and the B;'s are pairwise relatively 
prime. 
Since u, + vz; is odd, uy and vz must have opposite parity, i.e., uz, + Ve = 
1 (mod 2). This concludes the proof that (uz» Vz-) satisfies (1) for each ¢ € Tf. 
Finally, since S = 2u;(uz, + vi) = 2h,A, is independent of t, the k PPT's 
generated by (u,, v,), t¢ € T, have equal perimeters. 


3. Infinity of H, 


The main argument of this section rests on the following existential 
result; see [2], page 371. 


Theorem 3.1: For every positive number ec there exists a number &€ such that for 
each x, x > —, there is a prime number between x and (1 + €)x. (It will be used 
to prove the following proposition which has a certain interest in itself.) 


Proposition 3.2: Let k = 3 and let 6 > 0. Then there is a number € such that for 
every y, y > €, there are k consecutive primes Bi» By» re B, and a prime P, 
satisfying the inequalities 
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ys 8, < By <' +++ «Bye Vi + by, 


Poo) tea B,B 
Ganda E ca onesataec. 4 P, < (is BE Sie 2 ‘ 
By By 


scot Ties 


Proof: Let ¢, be a given number such that 0 < €, < v1 + 6 - 1. By Theorem 3.1, 
there is a number €, such that for every x > —,, there are at least k consecu- 
tive primes B,, Bo, ..., By in the open interval (x, (1 + €,)x). Let 
a 1+ 6 
(1 + €))2 


and take €) so large that for each ©, © > &5, there is at least one prime 
number in the interval (x, (1 + €)&x). 


Let € = max(&,, 2). Then for every y, y > —, we have that the interval 
(y, (1 + €))y) contains k consecutive primes, 


y < By < By < eee < B << (1 + €))y> (6) 
and the interval (y, (1 + €)y) contains a prime number Pax 

y <P, < (1 + edy. (7) 
We show next that the interval 


Ba Bin tsa, 1B BaP sire UE: 
X: Y] = Sthieileohinenicantt, (1 a 6th 
By By, 


contains P,. On_the one hand, we know from (7) that (xX, (1 + €)X]) contains’ at 
least the prime P,, since for k 2 3, X = ByB3 ... By/B, > Bz and Bp > y by (6). 
On the other hand, [X, (1 + €)X] is a subinterval of [X, Y] if we show (1 + €)X 
< Y. This last inequality is equivalent to 


B,B, jena By} Pe 56 BoB3 rate B, 


1+ 6 —E EE —————————— 
( ) B, (1 + €,)* By 





which, in turn, is equivalent to 
2pR2 2 
(i+ 6,)°8 2 > Bee 


But (1 + €,)B, > (1 + €,)y > B, by (6). Thus Y > (1 + €)X. This concludes the 
proof. 


We are now ready to prove the main proposition. 


Proposition 3.3: Let H,, k 2 3, be the set of integers S such that S is the per- 
imeter of exactly k PPT's. Then H, is infinite. 


Proof: Taking 6 = 1 in Proposition 3.2, we can count on k consecutive primes 
Pin Pas wees By such that 


B, < ¥2B,, 


so condition (2) is satisfied; moreover there is a prime P, such that condition 
(4) is satisfied. 
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Defining A,, u,, and v, as in Proposition 2.1, we see that (3) is also sat- 
isfied, so we may conclude that (u,, v,), t € T, generate kK PPT's having equal 
perimeter S = 2P,A,. 

Since y in (6) may be taken to be any number larger than €, it is clear 
that the above process may be repeated infinitely often. Each time we obtain a 
new set of k PPT's having equal perimeters. 

It remains to show that no PPT, other than the ones constructed, can have 
perimeter S = 2P,A4,. To do so, assume (u, v) generates a PPT with perimeter S$ 
= 2F,A,.- We will show that (u, v) is not a generator of a PPT unless (u, v) is 
one of the pairs (u;, v;) constructed above. 

Since S = 2u(u + v) = 2P,A, = 2B\Bp ... B,Py, there are but a finite number 
of possible values for u and u + v. We assume first that P, is a factor of u 
and consider the three possibilities: 


(i) u=P,, u+v = B By ... Bs 
(ii) u= BB, ... BP, utve=tl, 
(441) U = G1q0 +++ GmPes UF V = Gna1Ume2 9° Te 
where 47> +++ Gms me {1, 2, ..-, k - 1}, denotes any one of the products of m 
different primes from the set {B,, Bos sees By} and Im+1Um+2 °°* Uk the product 


of the remaining primes. 
In case (i), condition (4) implies 
2a QP, ADD, cae Beng lBy¢ BYBy 00s By H's 


so that u < v, a contradiction of (1). 
For case (ii), v = 1 - u < 0, contradicting (1). 


For case (iii), using (4), we write 


(4192 ++* Gm) (Gmai1Imee *** UIPy = AyP, > AZ/BF 


2 2 
BE nae BE 2 UirGige «oe Q)*: 


Then 
u= 9192 eee Vn > Ym+14Um+2 ate qx =u oe Cr 
contradicting (1). 


Next, we shall assume that P, is not a factor of u. Then P, must be a fac- 
tor of (u + v), and we consider the four possibilities: 


(I) utve=P,, u = BB eee Bus 
CEE) u+rv B,Bo eee Blas u™ Rs 


(TIT) U+ 0 = Quai Uma2 22° WPps U = W1U0 ©9* Um? 


where 414) +++ Ins ME {1, 2, ..., kK - 2}, denotes any one of the products of m 
different primes from the set {B,, Bp, ..., Bul, and 9.19.45 +++ Q, the prod- 
uct of the remaining primes. Note that x» + v contains at least two of the 
primes 2; as factors. 


(IV) ut+v = B,Pys u = ByBo «2+ ByuyBeq1 oe Bus t € T. 


In case (1), using (4), we get 
u+v =P, < 2A;,/B% < ByBp ... By = Us 


contradicting (1). 
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In case (II), v = B\Bo ... BP, - 1 > 1 = u, contradicting (1). 


For case (III), using (4), we have 


BP Gian +0 Gy > Syaitees yA ,/ Bt 
= (Yn+1 Um +2 eee I)? (4195 sae q,,) / Be > 2419 ore Vm = Zu 
a contradiction of (1). 


Case (IV) is seen to describe the k pairs (uz, v,) defined above. These k 
pairs then generate k PPT's with equal perimeters S = 2P,A,, and no other PPT 
can have this perimeter. 


4. Examples 


Let us conclude with a few examples. 


(1) When k = 3, we have: 


(165, 82) (143, 142) 81,510 
(1333, 628) | (1147, 1132) | 5,228,026 
(357, 118) (323, 202) 339,150 


(357, 194) (323, 286) 393,414 
(667, 158) (575, 382) | 1,110,550 
(713, 476) (667, 604) | 1,695,514 
(713, 534) (667, 666) | 1,778,222 
(1073, 198) (899, 618) | 2,727,566 








(2) Finally, let k = 4 and 
By = if. Bo = 19, B3 = 2s By = 23. 


For the integer P, within the bounds in (4), we can select any prime P, in 
the set 


{541, 547, 557, 563, 569, 571, 577, 587}; 
moreover, Proposition 2.1 allows us to take any nonprime P, in the set 


{545, 559, 565, 581, 583}. 
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ON A GENERALIZATION OF THE FIBONACCI SEQUENCE 
IN THE CASE OF THREE SEQUENCES 


Krassimir T. Atanassov 


Institute on Microsystems, Lenin Boul. 7 km., 1184-Sofia, Bulgaria 
(Submitted August 1985) 


A new direction for generalizing the Fibonacci sequence was introduced in 
[1], and [2]. In this paper, we shall continue that direction of research. 

Let C,, Co, ..., Cg be fixed real numbers. Using C,; to Cg, we shall con- 
struct new schemes which are of the Fibonacci type and are called 
3-F-sequences. Our analogy is of [1] and [2]; the form is 


Gg = Cys Dg = Cos Og = Cys a) = Cys Dy @ Ces oy = Ce 


= pl 1 
Fn+2 ~ n+1 * Yn (n 2 0), 


a: 2 
bisa * Baa, + Yh 
a 3 
Cn+2 Fh * Y, 


waste (ah ss vs #4) is any permutation of (@,,15 b,41> Cy4,) and (yl, y?2, 
y~) is any permutation of (a,, by, cy). 

The number of different schemes is obviously 36. 

In [3], the specific scheme 


Ay = Cy» by = Cos Cy = Cys a) = Cys dD, = C5, Cc) = Ce 


a = pb +e 
a + 
- 2 n+1 n (n > 0), 
n+2 ™ ns, + Fy 
Cn+2 % Qn] + b,, 


is discussed in detail. For the sake of brevity, we devise the following rep- 
resentation for this scheme: 


a@ b «e 
S= ‘ e \ (1) 
e @ »b 


Note that we have merely eliminated the subscripts and the equal and plus sym- 
bols so that our notation is similar to that used in representing a system of 
linear equations in matrix form. Using this notation, it is important to re- 
member that the elements in their first column are always in the same order 
while the elements in the other column can be permuted within that column. 
Every element a, b, and ec must be used in each column. 

We now define an operation called substitution over these 3-F-sequences and 
adopt the notation [p, q]S, where p, qe {a, kL, ec}, p = q. Applying the opera- 


tion to S merely interchanges all occurrences of p and q in each column. For 
example, using (1), we have 
e b a 
fa, cJS= {db a ep, (2) 
€ ¢@ 8 
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Note that in the result we do not maintain the order of the elements in the 


first column. To maintain this order we interchange the first and last rows of 
[2] to obtain . 


which corresponds to the scheme 


= ' = ' = ' = ' = ! = 
ay = Ci» by = Chs ey C3» ay Cy» by Ces Cy CE 
a =@¢e +p 
n+2 n+1 n (n > 0), 
n+2 Ane) + Cy 
Cn+2 ™ bn +) + a, 


where Cis Cos os CE are real numbers. 


We shall say that the two schemes S and S’ are equivalent under the opera- 
tion of substitution and denote this by 


5 + §', 


It is now obvious that for any two 3-F-sequences S and S', if [p, q]JS « S', 
then [p, g]S' «+ S. To investigate the concept of equivalence to a deeper ex- 
tent, it is necessary to list all 36 schemes: 


a «a a A “a *a G @ @ a a@ ea 
S,;=<4b b b So Db b @ S,=4b e¢ b Sy, *4b ¢ #2 
Cc Cc e@ ee Db ec be © iB 2D 
a “a 4 a qa WD ae Ia aD Ga @ 
Se"@Ab b @ ig=<> 2 «@ Sy "4b @ 2@ Sg*4b e e 
eo & e e «@ eb e bp a 

‘2 a e a e a ae ‘ a | 
S. =<) ib ia Si s=cp- & So ep © aS =<bp je @ 
ec eb sia GC ec a a © “bp "= @ ip a 
G@ tb “a a DB ai a ~b -@ a Db -@ 
S13 =<b a b Sin =<b a of Sis =<b @- Si6 =<b ec ec 
Cc Cc e Cc 2. 2D ec a e Ge 'e Dp 
a dD. Dd a Db 2D a Bb. & @ ~ iba 
S17 =<b a et Sig =<b a Cc Sig =e, GC 2 S56 = <b G¢ @¢ 
ec ee ce @G@ @ C @ e C-> 1 8 
a. Dee aD ze a «Db i a bp €e 
So, = Dp i@ “4 Soo =<b a b So3 =<>b @ a 5... =<b ec b 
e e b ey. .<@a © Ce @ Db e° 3@ <a 

we ie a e cf ‘ e st ' eC et 
Ss. =<p <a. B Sep ia) iC 5S. =m<p- DD b> S.. =eb DBD «@ 
3 C Be ” ed b i ec a Cc ” e -@..B 
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? e b ee b 

Sea <p a2 @ SS. = 25 a e s.. = 
29 

Cc b e ” e b a ” 


a 3€ e a ¢€¢ 6 
Syn “4b 2 a S3, =  @ Bb S35 = 
e 6b “db @ b @4 


Note that S = S53 and S’ = S39 so that S53 — Sa9- 


I 
—_——_—_—_ ——_— 
FRA QTR 
RFQ RTO 
rAd QgQRaeDT 
—/ WD 
Ww Ww 
a Np 
i) i 
es E_——_—_ 
QT 
RFQ 
go > 
—_—_—— 


QTR 
gro 
QR FQ 
—— 


(6 


We say that a 3-F-sequence S is trivial if at least one of the resulting 
sequences is a Fibonacci sequence. Otherwise, S is said to be an Essential 
Generalization of the Fibonacci sequence. 

Observe that there are ten trivial 3-F-sequences. They are 


Bis Sep Be Ba Ba 


a ee SG Ee 513» S17» 597» 536° 


These 10 schemes are easy to detect since they have at least one row all with 
the same letter. Furthermore, for these schemes one of the three possible sub- 
stitutions returns the scheme itself. For example, 


[b, elSi- Sis t+=1, 2, 3, 4 
las bISi< Saat by Seis 
[a, e]JSi + Si, t = 1, 10, 27, 36. 


The twenty-six remaining schemes are Essential Generalizations of the Fibonacci 


sequence. For eight of these schemes, the result is independent of the substi- 
tution made. That is, 


[p» q]S_ +* Sy [p> q]5\5  So5 
[p> q]5o9 + S33 [p, q]5.3 +* S39 


for all p and q. This means the substitution operation for these schemes is 
cyclic of length 2. 
For the other eighteen Essential Generalizations of the Fibonacci sequence 


schemes, all three possible substitutions generate three different schemes. 
For example, 


la, b]S, > Sz, la, b]Sg + S35 [a, b]Si¢ + S3y 
la, e]S, > Sy, la, e]Sg + Sy, [a, e]S,¢ + Sou 
[b, elS, + S,, [b, c]S, + 3), [bs c]S,6 + Soe 


All of the substitutions associated with the remaining eighteen schemes and 
their results are conveniently illustrated by the following three figures. 
That is, these pictures determine all possible cycles. 

For example, 


la, b]Syg +* Sjqg 
and 


[a, b]([b, el([as e]5,))  Soq- 
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FIGURE 1 FIGURE 2 





FIGURE 3 


Note that the figures tell us that many of the schemes are independent. 
That is, Sj, and S,g are independent. In fact, 5,g is related only to the six 
schemes listed in Figure 3. Similar results can be found for the other 
schemes. 

The closed form equation of the members for all three sequences of scheme 
S93 is given in [3]. By a method similar to that given in [3] or [1], the 
closed form equation of the members for all three sequences of the other 
schemes can be determined. We leave this task to the reader. Obviously, these 


results could be generalized to the case of four or more sequences with very 
Aiteie ditiiculty. 


References 
1. K. Atanassov, L. Atanassova, & D. Sasselov. "A New Perspective to the 
Generalization of the Fibonacci Sequence." Fibonacci Quarterly 23.1 
(1985) :21-28. 
2. K. Atanassov. "On a New Generalization of the Fibonacci Sequence." Fzibo- 
nacet Quarterly 24.4 (1986) :362-365. 
3. J. A. Lee & J. S. Lee. "Some Properties of the Generalization of the 


Fibonacci Sequence." Ftbonacct Quarterly 25.2 (1987):111-117. 


FREE 
10 


[Feb. 





ON THE PERIODS OF THE FIBONACCI SEQUENCE MODULO M 


Amos Ehrlich 


School of Education, Tel-Aviv University, Israel 
(Submitted May 1986) 


This article deals with Fibonacci's sequence 
ty * Op ty * Tete gy F my 
and with the arithmetical function 


K(m) = length of the period of Fibonacci's sequence 
when reduced modulo m. 


In the last few years I had some occasions to guide activities in a mathe- 
matics-with-computer club for 15-year-olds, where we investigated the function 
K(m). Theorems 1 and 2 of the present article were found (without proofs) by 
members of these clubs. To be more specific, these are those of the students' 
results, which I was not able to find in the literature either before or after 
they have emerged in the club. The rest of the students’ discoveries can be 
found either in [1] or in [4]. One of these is the following lemma which was 
suggested by the student Oded Farago. 


Lemma: For any m and n, K([m, n]) = [K(m), K(n)]. 


Proof: Follows from [5], Lemma 13. 


Theorem 2 in [1] says almost the same, but only for m and m that are rela- 
tively prime, so Oded's present version is more symmetric. (The lemma holds 


for every sequence that becomes periodical when reduced modulo a natural 
number. ) 


Theorem 1: For any fixed m let i,(n) = K(m"™*1)/K(m"). ‘Then: 


I. A,(”)|m for all n; 
Il. A,(”)[A,(% + 1) for all n; 


III. there exists ¢ such that A,,(n) =m for alln2t. 


This theorem emerged from the work of four girls: Shoshi Pashkes, Sigalit 
Teshuva, Mali Gana, and Chenit Lotan. 


Proof: 
(i) If p is prime and t is the largest integer such that K(p*) = K(p), then 
Theorem 5 in [1] implies 
Rte Fon <- t= 1 
Ap(n) = ’ 
pifn2t 


from which I, II, and III immediately follow. 
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(ii) If m = p®, then the conclusion follows from (i), since 
Am(™) = Ap(me)rA, (me + 1) ... Ap(me +e - 1). 


(iii) Now let (a, b)=1 and assume that the theorem holds for m = a and m = b. 
By hypothesis and the lemma, 


hg(m)|a, 4p (m) |b, K(a"b”") = [K(a"), K(b")], 
also 
K(a"*1p"tl) = [\,(n)K(a”), A, (n)K(b")]. 


Let \g,(") = K(a"t1b"*1)/K(a"b”). Then A2,(2)|Ag(“)Az(m). Thus, AQg,(") |ab. 


Let p be a prime such that p*” 
ity, let pia. 


Ag(n), p**||K(a”"), p?||K(d). 
Since pb, by part I we have p/A,(m), so p°|K(b”) for all n. Therefore, 


2, = Maxix, + y,,. c} - Max{y,», c}> 








Agp(m). Then p|ab; without loss of general- 





Let p*"| 





that is, z, = 2,» X, t+ Y¥, ~ Cs or 0. By hypothesis, 2x, < 2,4), and Y, < Y¥n41- 
Therefore, 2, < 2,41» SO p*"|Agp(m + 1). Since p is arbitrary, we have 


hap(n) |Aan(m + 1). 


By hypothesis, there exists tq such that Ag(”) =a for all 7 2 tg. Since 





Aqg(") = a means that K(a"*!) = @aK(a”") this implies that y,4) 2 y, + 1. It fol- 
lows that there exists a T > t, such that for all n 2 T we have y, > c and thus 


S, = 2y- 


Since, for such an 7, \,(n) = a, it follows from 2, = x, that p*"|la. 
Since p{b, it follows that, for all n => 7, p**|lab. 


For ” sufficiently large, this holds for every prime p that divides ab; for 
such an 1, hgp(n) = ab. 


Theorem 2: For any even 7 >3, K(u;) = 2¢. For any odd z > 4, K(u;) = 4. 


Remark 1: Amihai and Moshe, the boys who found this, used different words. 
They said that the elements of the sequence K(u,)> K(us)> K(ug)> ‘ee GEG, 
alternatively, the elements of two arithmetical sequences, one with the dif- 
ference 4 and one with the difference 8. 


Remark 2: The second part of Theorem 2 follows from Theorem 3 in [3]. 


Proof: K(m) is the first 7 after 0 such that u; = 0 and u;+) = 1 (mod m). 

Theorem 3 in [1] says: For every m there is ad such that u; = 0 (mod m) 
if and only if d|j. 

If m =u; > 1 then d = t, since the elements before u; are not changed when 
reduced modulo u;. [This proves that K(u;) is a multiple of 7.] 

For the same reason, if ~ > 3 then u;-, # 1 (mod u;). 


Now uj+] = uz-] (mod u;); therefore, if 7 > 3, then the i‘ element of the 


Fibonacci sequence modulo u; does not start a new period, instead, it starts a 
sequence of u;-,; multiples (mod u;) of the original sequence. Hence, 


=— = 2 
Ugi41 = Ugz-1 = UZ_, (mod u;). 


For every 7, us 1 = Uz_pu, + (-1)*; therefore, if z is even, then Ure41 = lL 
(mod u;) and K(u;) = 27. 
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For odd 7, uUy;,,; = -1; therefore, Unig] = 1, SO U3;,, # 1 and K(u;) = 4. 
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HOW MANY 1'S ARE NEEDED? 
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13508 Bartlett St., Rockville, MD 20853 
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Let f(”) denote the number of 1's necessary to express 7, using the opera- 
tions + and x (and parentheses). Determining f(”) is an old problem, origin- 
ally considered by Mahler and Popken in 1953 [3]. We have calculated f(n) for 
n < 310 and we present some statistics. (The reason for a power of 3 will be 
explained later.) 


The Problem 


With f(”) defined as above, we have 


( l, for n= 1, and 


fim) = (1) 


ab=no 
a+b=n 


min {f(a) + f(b)}, for n> 1. 


This formula is very time consuming to use for large n, but we know of no 
other way to calculate f(n). 

The behavior of f(n) is interesting. Selfridge has shown that 3% + Q3k-1 
is the largest n for which f(n) = 3k + 0, for 0 = 0, +1. His proof is by 
induction, and the induction step is based on the following observation: If 
b(m) is the largest m for which f(n) =m, then b(m) is the largest element of 
the set 


U {b(r) + b(s), b(r)b(s)}. (2) 


r+ts=m 


Using (2), it is fairly easy to show that b(m) has the required form. 
There are two competing conjectures about the behavior of f(m) for large n. 
It has been conjectured that 


f(t) < 30. + e)log,n, for large n, and any e > 0. (3) 


It has also been conjectured that there is a set S (possibly of positive den- 
sity) and a positive constant ¢ so that 


f(n) > 3(1 + e)log,n, for all n in S. (4) 


The Results 


We calculated f(n), using equation (1), forn < 310 (59,049). We chose 310 
because we would have all n for which f(n) < 30, by Selfridge's results. 

We broke the interval into 30 subintervals between the values of the form 
3k + 93k-1 for 6 = 0, +1 and we also looked at the sets S(m) = {n: f(n) =m}, 
for m= 1, 2, ..-, 37, incomplete beyond m = 30. 
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Analysis of the thirty subintervals 


One typical subinterval is the interval 18 < n < 27. The values at the 
"endpoints" differ by 1. If conjecture (3) were true, we would expect the val- 
ues of f(n)/log.n in the interval to "flatten out" and approach the values at 
the endpoints, as nm gets large. Table 1 gives the mean and standard deviation 
of f(n)/1og 3 in each interval. 

While the analysis of such small values of f(n) has very little to do with 
behavior at large m, it is clear that in this range conjecture (4) is strongly 
supported. 


The single worst value of f(m)/log,n encountered was at m = 1439, with f(n) 
= 26, and f(n)log,n = 3.9281. 


TABLE 1 
Mean and Standard Deviation of f™)/log3n, fora<n<b 





Analysis of the sets S(m) 


Let S(m) = {n:f(n) = m}. In Table 2 we consider the following questions 
about S(m). 


What is its first element? 
How many elements are in S(m)? 
What is its last element? 

What is its average element? 


One result about the sets S(m) not captured in Table 2 is: If b(m) and b,(m) 
are the largest and second largest elements of S(m), then b,(m) = [(8/9)b(m)], 
wheie [+] is the greatest integer function. 


Outline of proof: The proof is by induction. 


The result is true by inspection for small values of m. For large values 


of m we have an equation similar to (2): b,(m) is the second-largest member of 
the set 
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U {b(r) + b(s), b(r)b(s). b(r) + by (8s), b(r)b,(s)}. (5) 


r+s=m 


For m > 9, it is easy to show that both b(m) and (8/9)b(m) belong to this set. 
It only remains to show that there are no elements between these two values, 
and this can be done by a simple case-by-case examination using the results of 
Selfridge and the induction hypothesis. 


TABLE 2 
Analysis of the sets S(m) 


108 
162 
243 
324 
486 


729 
972 
1458 
2187 
2916 
4374 
6561 
8748 
13122 
19683 
26244 
39366 : 12534.0 
59049 17315.0 
78732 
118098 
177147 
236196 
354294 
531441 
708588 





Comments 


In his paper [2], Guy relays some questions about the function f(m). We 
comment on three of these: 


Q@: For what values a and b does vi"9") = 2a + 3b? 
A: (273°) = 2a + 3b for all 273% < 319, at least. 


OC: tt £(2#°s") = 2a + 3b and there is a larger n' so that f(n') = 2a + 3b 

(a, b > 0), must n’ = 273", for some r, s? 
A: No. Two counterexamples are 27 < 395, with f(27) = (335) = 14, 

and 2732 < 355, with f(2732) = f(3°5) = 20. 


Q: When the value of f(n) is of the form f(a) = f(b), with a + b =n, and 
this minimum is not achieved as a product, is either a or b equal to 1? 
A: Yes, at least for n < 3/9, 


The calculation of f(m) was performed on a Symbolics 3645 LISP machine using 
equation (1), and we used over 50 hours of CPU time. 
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1. Introduction 


Consider WV points placed on the circle of unit circumference in the follow- 
ing way: begin by placing a point anywhere on the circle. Now place another 
point so that the angle (or circumferential distance) between the two points, 
measured clockwise from the first point, is equal to a. The third point is now 
placed at a clockwise angle of a from the second point. Thus, we successively 
place W points on the circle by our angle a. 

Our problem is to find the value for aso that these points are _ spread 
about the circle in the most even (which we call optimal) fashion. We show 
that, in certain senses, the golden section (a = T = (5 - 1)/2) provides the 
optimal spacing of points, where the number of points can assume any value. 

This problem originally arose while investigating the phenomenon of phyllo- 
taxis—regular leaf arrangement. Most higher-order plants exhibit a remarkable 
degree of regularity in the positioning of their leaves. In a sunflower, for 
instance, one can perceive two sets of opposed spirals which each partition the 
set of florets. Intriguingly, the number of spirals are almost certainly 
consecutive members of the Fibonacci sequence 


y= Kolo eh Lynden Se 2yehyj = 20s, = 


This pattern (which we call Fibonacci phyllotaxis) manifests itself in 95Z of 
those plants which produce their leaves sequentially. In parallel to this 
observation, the divergence angle subtended by consecutively formed leaves is 
quite close in value to the ratio of these consecutive Fibonacci numbers. In 
the limit, F,,_,/F, is equal to the golden section. To simplify the situation, 
we consider just the angular displacement of the leaves and thus we develop a 
simplistic model of plant growth with leaves appearing as points on a 
meristematic ring, successively placed at a constant angle. 

What this paper shows is that the plant places its leaves in the optimal 
manner—in order to spread its leaves most evenly (and thus reduce leaf 
overlap) the optimal divergence angle is shown to be the golden section. The 
partition of the circle by the golden section is also examined in detail to 
reveal a rather self-similar structure. 

We use results from The Three Gap Theorem (originally the Steinhaus Conjec- 
ture) which states that the above WV points partition the circle into arcs, or 
gaps, of at most three and at least two different lengths! The result is all 
the more remarkable since it holds for all irrational a and for any number of 
points. It also holds for rational a = p/q with the number of points less than 
q- (For W = q the circle is partitioned into g equal gaps.) Even though this 
has been proved by various mathematicians ([1], [2]-[7]), the result does not 
appear to be well known. 

Note that, in order to conserve space, where complete proofs of results are 
not presented we either refer the reader to an existing proof or briefly 
outline a proof. 
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2. The Three Gap Theorem 


Suppose that we have consecutively placed NV points on a circle by the angle 
a. Let (u,(V), ug(V), ..., uy(NV)) be the sequence of points as they appear on 
the circle, ordered clockwise from the origin u,(NV) = 0. That is, 

{u, (WV), u2(V), ..., uy(W)} = {0, 1, 2, ..., MW - 1} where {uja} < {uj.,a}. 


Thus, for example, with a = ¥2 the first 12 points placed on the circle appear 
in the order (0, 5, 10, 3, 8, 1, 6, 11, 4, 9, 2, 7). We call u34)(NV) = ujey 
the successor to u;, or uj+, = Suc(u;). Equivalent to the original statement 
of The Three Gap Theorem is the fact that the difference between succeeding 
points assumes at most three, and at least two, different values. The 
following determines the ordering of points around the circle. (For a proof, 
see van Ravenstein [7, Theorem 2.2].). 


Theorem 1: Us O<sm<WN- uu; 
Suc(m) - m= ug -Uy, N-ug sm< Un» 

“Uy> uy <m< WN. 
Thus, for our example with a = V2 and WV = 32; 


x Olm< 7; 


Suc(m) - m= 


eo 2 Sines 22. 


It is easily seen that u; = 5(j7 - 1) mod 12, where y mod x = y - x[y/x] = 
x{y/x}. In general, if NV = uy + uy, the circle is partitioned into gaps of 
just two different lengths and then 


uz = (C7 - 1l)ug) mod WN, jf = 1, 2, coos MN. (1) 


It is easy to see that the length of the gap formed by point m and Suc(m) 
is equal to {(Suc(m) - m)a} where {x} denotes the fractional part of x such 
that x = [x] + {x} where [x] is the largest integer not greater than x. In 


fact, for gap lengths less than 5 this gap length is equal to 


| (Suc(m) - mall, where |x| = min({x}, 1 - {x}) = |x - [x +$]|, 


the difference between x and its nearest integer. (This is always the case for 
N > q, [notation defined in Theorem 2]; in what follows, we will always make 
this assumption. Note that g, is the first point to replace 1 as the closest 
point to the origin.) Thus, Theorem 1 shows that the circle of W points is 
partitioned into N - up gaps of length |u zal, V - uy gaps of length |u,al and 
Uz + uy - N gaps of length ||u2al| + ||wyall. The same applies for rational a, say 
a= p/q in lowest terms, where WN < q. In this paper, however, we will always 
assume that q is irrational. 

Point uz is the successor to 0, while 0 is the successor to uy; that is, up 
and uy, are the points which neighbor the origin. We see that we need only know 
the values of these two points to determine the entire ordering. 

We can characterize the angle a by the following. Let V(x) denote the path 
of a defined to be a sequence of pairs (u2, uy), the points which neighbor the 
origin as points are successively included on the circle. For example, 


N42) = CO, Di. Gs Bh. Bs D, Ge-Jie Ch 10a). 


It can be shown ([7, Proposition 4.2]) that each point always enters one of 
the larger gaps. Two gaps are formed, one equal in length to the smallest gap 
present. Thus, it is natural to define the ratio of gap division as the ratio 
of the smallest to the largest gap present. Hence, we let 
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r, (a) = 
VO Tusal + Fava 
[This is in fact the ratio point W - 1 that divides some (large) gap.] 

The path sequence /(a) and the ratio of gap division r,;(a) are quantities 
we will use in our analysis of the golden section's unique distribution 
properties. We can in fact determine explicitly their values in terms of the 
continued fraction expansion of a, which is expressed by 


1 
a= G9 °+ 


a, + 
a, + 


1 
az eee 
{293 A,» Ans Ags «ee. 
The n*® tail of a is 
$, © (ays: Qay0 Syaae ose) (2) 
such that 
a = {ao A)» Ans ses Ay). tate 
We say that a is equivalent to B if some tail in a is equal to some tail in 8. 
Partial convergents are defined by the (irreducible) fractions 
Pnii _ Pn-2 + tPn-1 


+e {a a a eee a t} c = 1 2 eee a 
° 0? 1? 2? > n-l? > > ’ ’ n? 
in,t dn -2 t4n-1 





where 
Pn, an Py, 
——————— 2 9 = = 0; = = ‘ 
Pick: Tn Bao?" Gay ig Pay 1 


We call p,/q, a total convergent to a. 
The reader is referred to [7, Theorem 3.3] for a proof of the following. 


Uu 


Theorem 2: 
; 


Yn-19 n odd, <3 n odd, 
Uy = 


q 


s n even n even 
ase <i? > , > 


where q,,;-, << i < &, & 22). 


For qn-,; < N < qy,, (m 2 2); 


a n odd, . n odd, 


Un,i? 


In-22 m even, Qn-12 ” even. 


For W < q), uj = Fo =le tech ee scieus 

The following proposition may be easily proved from the definition of ry (a), 
Theorem 2, and the continued fraction theory. 
Proposition 3: 1 
i'+ t,’ Qn-1 <M 5 4n,15 


1 
Sota tent FP 


ry (a) = 


20 
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where 2 = 2, 3, ...5 Gy, (nN 2 2). [t, is defined by (2).] 


From Theorem 2, 


V(a) = (C1; q\)'> (Ques a t 


'; 23 eces And nN = 25 3s Saana (3) 
where 


(1, q1)> O0<a<35, 
(1, q,)' 


(q,> 1); s<ac<l, 


; ‘a agen n odd, 
Gn,i? Wn -] a 


(4,39 In-j)» m even. 
3. The Golden Section 


For convenience, let the partition of the circle of unit circumference by 
the successive placement of points 0, 1, 2, ... by the golden section, T, be 
denoted by G. The partition by 1 - t we denote by G’. 

The continued fraction of t is given by 


+= (& 1.4 t)-@ (66-3, 2: +} @ 406 1,.14. 153. 
All convergents to T are total convengents and 
Fe dgng Pig @ Egcg * Fg MW @ li Fy © 1a 2g Ds 


That is, convergents to t are equal to the ratio of consecutive Fibonacci num- 
bers. From Theorems 1 and 2, for F, < WV < Fy4j;5 


FPF.» OsmsN-FP,;> 


| Naot me < Fy ae n odd, 
“Foy F,.y Sm< Ws 


Suc(m) - m= (4) 
Fi,uzs Os ms Fm Fas 
OF eg eR ca es n even. 
-F > PS mee, 


When VV = F.4,» from (1), 
uj = ((-1)"-!(j - 1)F,) mod F,,,- 


Since F,,-, - F,,1 =(-t)” (by induction), (4) shows that W points (F, < N < Fy4}) 
partition the circle into W - F,-_) gaps of length al, N- F, gaps of length t” 
and F,+ 1 - N gaps of length 1”~4 


From Proposition 3, sincet,=1+1T, n=1, 2, 
Wit} = te b= BE «8 s F,,,. (5) 


Theorem 4.1 from [7] describes the partition G by looking at the transfor- 
mation of gap types as points are included on the circle. Gap types are either 
"large" or "small" when N = up + uy, that is, when V is the denominator of a 
convergent to a. For the golden section, this theorem describes the following: 
each large gap present when WV = F, is divided by the addition of a further F,,_, 
points into two new gaps which can be labelled (in clockwise order) as 
small:large (m odd) or large:small (m even) when W = F,4).- Those small gaps 
present when WV = F, can then be labelled as large. 
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This fact and (5) can be used to prove the following self-similarity prop- 
erty of G, thus demonstrating the beautiful symmetry inherent in Fibonacci 
phyllotaxis. 


Theorem 4: Consider any large gap present on the circle partitioned by the 
placement of F, points by the golden section. Include further points on the 
circle and observe the resulting partitioning of this gap. If we pretend the 
gap is itself a circle of unit circumference (by lengthening it by the factor 
t"-2 and identifying its endpoints as the same) then its partition is identical 
to G if n is odd, or equal to G’ if n is even. 


Let us interpret WV to be a time variable and define the age of a gap to be 
the time it has survived without being divided. That is, the age of the gap 
with endpoints u;, uj+; is V- 1 - max(uj, Uj+1)- From [7, Proposition 4.2] 
each point, for all a, divides the oldest of the larger gaps. Using [7, 
Theorem 4.1] it can be shown that only for the golden section does the 
formation of a large gap always coincide with that of a small gap. This proves 
the following. (Note that we assume that a is between 0 andl. If a >1, the 
following results hold if a is replaced by its fractional part, {a}.) 


Theorem 5: For 5 <a< 1, each point always enters the oldest gap if and only 
ifa-=t. For 0 <a< 45, each point enters the oldest gap if and only if a = 


t2 = (3 - ¥5)/2. 


_ 


Intuitively, in terms of phyllotaxis, it seems sensible that points be in- 
serted in the oldest gap as the above result shows. This property must ensure 
an ideal distributtion of points. In fact, the following theorem shows that the 
golden section provides the optimal value for gap division (our criteria for an 
optimal distribution) in the sense that the smallest value assumed by the ratio 
of gap division is largest for the golden section (where 5 <a < 1). However, 
the golden section is somewhat of a compromise as Theorem 7a shows (that the 
ratio of gap division's maximum value is smallest for the golden section, where 


I 2 
a % = 
oe eg). 


Theorem 6: , max min Yr, (a) = 12, exclusively attained by a =T. 
Swed 
max | min Py (a) = 12, exclusively attained by a = 2, 


O<a<> WN 
Proof: We first consider the case where ;< a < 1. From Proposition 3, 
1 i 


min Yr,(a) = min min r_(a) = nin —— = ——_, 
N . N Qn-1<NS4n N aioe i + max bn 


where n = 2, 3, ..- (q, = a, = 1 since a > +). 


Consider a =a’ #1, which has a, >1 for some integer k greater than l. 


Then max t, > t, > 2, sor (a') <1. The result follows since Lepr (rt) = 12. 
ve n k N 3 3 N 


The second statement follows by symmetry (note that r(1 - a) = ry (a)). 


Theorems 7a and 7b follow from Proposition 3 in a similar fashion. 


Theorem 7a: , min, max Pry (a) = 2, exclusively attained by a =T. 
27°a< 5 9 
min | max Py (a) = 14, exclusively attained by a = 2, 
N 
2 


Theorem 7b: , min max r,(a) = +2, exclusively attained by 
7s a< 1 N 


+ 


i: 


ate, iy¢. b+ 1} ==> 


22 


= 7? where qa is any integer greater than l. 


Q 
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min max Py (a) = 1%, exclusively attained by 
O0<a<4 N 
1 


a={0;a+1, 1+t}= 277 where a is any integer greater than l. 


We determine the value of « which ensures a path which consistently maxi- 
mizes the length of the smallest gap on the circle. For each value of WN, the 
points are generated by a constant angle a. This value of a may change with WV 
but only in such a way that the path is retained: so that the addition of extra 
points does not alter the relative order of existing points. We show that V(t) 
is the path which ensures that the smallest gaps are consistently as large as 
possible [where, initially, (u2(4), u3(4)) = (2, 1), 3 <ar< <. That is, as WV 
increases, if the pair (uj(V), uy(V)) does not assume the value equal to the 


appropriate successive element of V(t), then the smallest gap thus formed will 
not be as large. 


Note that the golden section has path 


V(t) = (li, 1D); (2s Ts C25. Ss C55. 3) 5 oue)is 


” (Cl, B) Peet? Fars a = 24 34 ere 
Theorem 8: Suppose that (u2(4), u3(4)) = (2, 
a where >< a <5- Then V(t) 
length of the smallest gap. 


1) generated by a constant angle 
is the path which consistently maximizes the 


Proof: We prove the result by induction. Initially, ; <a< 3 or a= {0; 1, l, 
t3}, 1 < t3 < , such that point 2 is closest to the origin. The next element 
in the path must, from (3), be (2, 3). From Proposition 3, point 2 is furthest 
from the origin if a3 = 1 than if a3 > 1 since then it divides the gap bordered 
by the origin and the first point into a larger ratio. Hence, a = {0; l, l, l, 
ty}. This ensures, from (3), that the next element inthe path is (5, 3). Thus, 
the first three terms in the path belong to V(t). 

Now, assume that the terms in the path equal successive Fibonacci pairs and 
that (u2(N), uy(V)) = (Fy-15 Fn) where F, < NV < Fy4,], n even. Then, a = {0; 1, 
1, 1, ..., 1, t,} (n - 1) ones, 1 < t, < ~. The next element in the path must 
be (F,,41, F,) succeeded by (F,4 1,5 Fy42) if a, = 1. From Proposition 3, the small 
gap bordered by origin and point F,, is larger if a, = 1 than if a, > 1. The case 
is similar for odd nm. Thus, the path is equal to V(t). 


Note that the theorem shows that maximizing the length of the smallest gap 
ensures convergence to the golden section. Similarly, V(t?) consistently maxi- 
mizes the length of the smallest gap where, initially, there are three points 


on the circle and i< a <3. The following generalizes Theorem 8. Its proof 


is similar in manner and is omitted. 


Theorem 9: Suppose that we have placed q, + 1 points (m 22) generated by a = 
{0; a5 Go. «++» A,» tyy,}. Then as more points are added, V(a’), where 

Pn * TPy-1 

In + TI-1 

is the path which consistently maximizes the length of the smallest gap. 


a’ = (03 ays Qos econ Qys 1+ 1} = 
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1. Introduction 


Let k be a positive integer and let {T7,}°_, be a kth-order integral linear 
recurrence defined by 


Path = MWnapey + Folngp-2 + °° + GF, (1) 


with arbitrary initial terms 79, T,, ..., Ty-,- Associated with the recursion 
relation (1) is the characteristic polynomial 


f(x) = xk - aa"? ee ie (2) 


with characteristic roots Y}, fo, ..., ry. We will seek subsequences of {Ty} 
such that the recursion relation (1) is also satisfied as a congruence modulo 
some integer m. Specifically, we will endeavor to find positive integers d and 
nm such that 

Prekd = FTn+(e-1a + F2Tna(k-2.a Fett M-1Tn+g + UT, (mod m) = (3) 


for all nonnegative integers nm. This investigation was suggested by Freitag 
[2] and by Freitag and Phillips [3] and [4], and will generalize the results of 
these papers. 


Two approaches will be taken in satisfying congruence (3). In the first 
approach, given a fixed modulus m we will seek to find integers d such that (3) 
is satisfied. Along these lines, Freitag [2] proved the following theorem: 


Theorem 1: Let {F,} as usual denote the Fibonacci sequence. Then 

Fisoa = Faca + F, (mod 10) (4) 
for all nonnegative integers n if and only if d = 1 or 5 (mod 12).0 

The second approach will be to take the integer d from among the integers 

appearing in a specified sequence such as the sequence of primes and then find 
moduli m, depending on d, such that congruence (3) is satisfied. Corresponding 
to this approach, Freitag and Phillips proved Theorems 2 and 3 in [3] and [4], 
respectively. 
Theorem 2: Let {T7,,} be a second-order recurrence defined by 

T+2 = 242 q41 + AT. 
Then, if p is a prime greater than 3, 
Tn +2p = n+p + apT,, (mod 2p) 
for all nonnegative integers n.{J 


a,T 


Theorem 3: Let {7,} be a k*®-order recurrence with distinct characteristic roots 
satisfying 


Suet = QT, 4K-1 + Aol 4K-2 + vee + a,f,, « 


Then, if p is a prime, 
1989] 25 





CONGRUENCE RELATIONS FOR k*"-oRDER LINEAR RECURRENCES 


Trekp = 17 ne(k-vp + FTne(e-2)p ttt + G-1Tntp + %Tn (mod p) 
for all nonnegative integers n.{J 


2. Definitions and Known Results 
We will need the following definitions and lemmas to continue. 


Lemma 1: Let {7,,} be a k""-order linear recurrence with distinct characteris- 
tic roots r), Po5 +++» Y,- Then 


m 
T, = ae (e{? + cn + eee + aye Rtas 


t t 
t=1 


where the c are complex constants and s; is the multiplicity of the root r;. 


Proof: This is a classical result in the theory of finite differences (see, for 
example, Milne-Thomson [5, Ch. XIII]).0 


Definition 1: The primary linear recurrence Gil ce is the recurrence satisfy- 
ing (1) and defined by 


a Saiee n n 
VS rp tra tees + ry, 
where Y), Yo, «++, Y, are the zeros of the characteristic polynomial C2) EF 
any characteristic root 7 * 0, we define r? to ibe i. 
Lemma 2: Suppose {7,,} is a k*®-order linear recurrence satisfying 

Taek = UTInek-1 + TTygp-2 Free + ATy- 
Suppose m is a positive integer such that (a,, m) = 1. Then {7,} is purely 
periodic modulo m. 


Proof: This is proved by Carmichael [1, p. 344].0 


Lemma 3: Let {T7,,} be a k*®-order integral linear recurrence with character- 
istic roots Py> Poo veer Pye Let h be a fixed positive integer, and let q be a 
fixed nonnegative integer. Then the sequence 


{Suhpeo = (Tinsqhn=o 


also satisfies a linear integral recursion relation 


_ oh) (A) ¢ (n) 
Snek = 7 Snan-1 + 22 Spag-2o toe + MR Sans (5) 


where ar, as”, pepats ay” are integral constants dependent on h but not on q- 


Further, if gy is a fixed integer such that 1 < jg < k, then 
(h) _ eras Sits h 
a; > (-1) Py Py cee Me (6) 


where one sums over all indices t)> tos erate 1; such that 


1st, <4, < +++ < tj sk. 
Proof: This is proved in [6].0 
3. Main Results 


We now present our principal theorems. 


26 [Feb. 





CONGRUENCE RELATIONS FOR kt -ORDER LINEAR RECURRENCES 


Theorem 4: Let {T7,,} be a k‘"-order recurrence defined by 
Prek = UT nay + Tra x-g + oe* + AyTh> 
Let p be a prime. Then for all nonnegative integers b, 


Ln+kp> = Ql +(k- 1) p + Aol + (K-2)p? + eee + Ay 12 4p! + a,T,, (mod p), 


where ” is any nonnegative integer. 


Proof: Let Ps Vos +++s Ym be the distinct characteristic roots of AT }s Let 
denote the integers of the algebraic number field Q(r,, rj, .--, %,), where @ 
denotes the rational numbers. Let Z denote the rational integers. Let P be a 
prime ideal of R dividing p. Let o be the Frobenius automorphism of the finite 
field R/P having Z/p as a fixed field. Then o is defined by o(x) = x?. Then, 
for any nonnegative integer b, 0%, defined by o9(x) = x”, is also an auto- 
morphism of R/P fixing Z/p. 
Now, for i <4 <im; 
ne = “a? + a,er"* tees tar, + a. (7) 


Applying o” to equation (7), we have, for 1 <7t<™m, 
kp? . b pk 
a*(s*) = rxP = o*(a,7t) % a3 + see + G,) = pes? (P; 
k . 
=: } a, rfk-5)P? (mod P). (8) 
j=1 


By (8), (1), and Lemma 1, we have 


m 
pa = b 
Tecate * ZL [(e + Dn + eee + on™ ')r? |e}? 


t=1 


t t 


m By - k = b 
’ ze (O) 4 en ius a ‘et hag nik Dp 


t 


: st 3 


m . , 
a; Ue 00) + Un +... + gir gan 


k 
= 22 a; ar, jp? (mod P). (9) 


Since the first and last terms of (9) are rational integers, we have 


k 
Pa+kp> = 2 aj T na ck - 50? (mod p).U 
j= 


Remark: We note that Theorem 4 is a generalization of Theorem 3. 


Theorem 5: Let {Ty} be a k*®-order recurrence defined by 
Taek = Taney + Fotngx- ae et aT,» 


Let ¢ be a fixed positive integer such that (c, a@,) = 1. Then there exists a 
fixed modulus g such that if h = 1 (mod g), then 


Taekn = TT nece-rn + F2Inece-2n ttt + MH-Tnen + UT (mod e), 


where ” is any nonnegative integer. 
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Proof: If h is any positive integer, then by (5) and (6), 


=» A) (h) (h) 
Treen = 21 Inace-vn + 42 Tnoce-2yn t °° + KR Tne — 
where, for 1 <j <k, 
eA Lh h 
y> (-1) WP cee Tye (11) 


where one sums over all indices ty» to» deen Be 


4 such that 


lst, <t, < o++ < tj sk. 


Let n; = (5) . Let 1 < j < k be a fixed integer and let ; i". areiers i de- 

note the (‘ ) algebraic integers (-1)7*1y, Yr, «+.Y, 5 where these represent all 
. 1 2 dj 

the (*) products taken j at a time of the characteristic roots Y), Yo. «++» Yr, 


of {7,}. By the theory of symmetric polynomials, for a — integer j such that 


wr hag sees i are the roots, possibly 


1 < j<k, the n; algebraic integers t 
with repetitions, of a monic polynomial of degree n; itt rational integral co- 


efficients. 
Let ‘v4. defined by 
v= (t -5" oe ay" sales <r 


be the primary linear recurrence with characteristic roots “? : ‘, eiats ef. 

Since (a,, ¢) =1, it follows by Lemma 2 that {V, ) is purely periodic steko c. 
Let d; denote the period modulo ec of {vy } for 1 <g<k. Let g be the least 
common multiple of d,, dj, «...s dx. Since by (11), 


ys a) an a. : 


d 


it follows that if h = 1 (mod g), then 


ay = . * = y = a; (mod c). (12) 


The result now follows by (10).0 


Corollary: Let {7,} be a k**-order linear recurrence defined by 


45 '* Ql, 4k- 1 Gola xeg t+ *** + ag%,. 


Let p be a fixed prime such that p] a,. Then there exists a fixed modulus g 
such that if h = pb (mod g), where b is any nonnegative integer, then 


Tekh = UTnece-pn + F2Trace-2yn + °°* + %2_ (mod p), 


where m is any nonnegative integer. 


Proof: Let {V,} be any primary linear recurrence witn characteristic roots Yr), 


Po: ea > Pr Then 
oar p> p> . Ps p> = 
ie: ik Oe eS eee a | 
(mod p). 


Let the primary linear recurrences {v} and the integers a, wherel1<sjg<k, 
be defined as in the proof of Theorem 5. Choose the modulus g in the same man- 
ner as in the proof of Theorem 5, letting p=ec. Then 
(9) x: 9) 
Voo = V, (mod g) 


and 


a™ = a = VD = = vy = 


j a; (mod p) 
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for all j such that 1 < jg < k. The proof now follows by (10).0 


Remark 1: Note that if p is a fixed prime, the corollary to Theorem 5 is a 
strengthening of Theorem 4. 


Remark 2: Theorem 1 follows from the corollary to Theorem 5. By the proof of 
this corollary, it can be shown that if d = 1 or 5 (mod 12), then 

Faaa = Fag ® f, (mod 5). (13) 
Similarly, it can be shown that if d = 1 or 2 (mod 3), then 


Pisoa = Fn+oa + F, (mod 2). (14) 


It thus follows that if d = 1 or 5 (mod 12), then (14) holds. Since 2 and 5 
are relatively prime, it follows from (13)-(14) that if d =1 or 5 (mod 12), 
then congruence (4) holds. This proves the necessity of Theorem 1. The 
sufficiency of Theorem 1 follows from the fact that {F,} has a period modulo 10 


equal to 60. Examining (4) for all integral values of d between 1 and 60 es- 
tablishes the result. 


Theorem 6: Let {7,} be a k*®-order linear recurrence defined by 
Trek = UWI g-1 + UThee-2 tore + &T,- 


Let ¢c be a fixed positive integer such that (¢c, a,) = l. Then for all non- 
negative integers b, there exists an infinite number of primes p of positive 
density in the set of primes such that 


Prokp> = Tne (n-1ype + F2Ins (e-2yp> Foe? + S-1T ag pot UHTn 


(mod ep), (15) 


where ” is any nonnegative integer. Furthermore, there exists a fixed modulus 
g such that if p = 1 (mod g), then congruence (15) is satisfied. 


Proof: By Theorem 4, the congruence (15) is satisfied modulo p for any prime p. 
Given the integer c, we choose the modulus g in the same manner as in the proof 
of Theorem 5. By Dirichlet's theorem on the infinitude of primes in arithmetic 
progressions, there exists an infinite number of primes p such that p =1 (mod 
g). Further, the density of such primes is 1/$(g), where $ denotes Euler's to- 
tient function. By Theorem 5, congruence (15) is also satisfied modulo e, 
since p> is also congruent to 1 modulo g for any nonnegative integer b. Since 
we can also assume that (p,c) = 1, it follows that (15) is satisfied modulo 


ep.U 
Corollary 1: Let {7,} be a k**-order linear recurrence defined by 

Letk * Sinek<t © Seigagae * *** t Ugly 
Let @ be a fixed prime such that cla,. Then for all nonnegative integers b, 
there exists an infinite number of primes p of positive density in the set of 
primes such that 

Thekps = Tne ce-ve® + F2Tne(e-2pe tote + Ay Tnepe + Ul 

(mod ep), (16) 


where m is any nonnegative integer. Furthermore, there exists a fixed modulus 
g such that if the prime p = c> (mod g), where b is any nonnegative integer, 
then congruence (16) is satisfied. 
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Proof: This follows by the corollary to Theorem 5 and the proof of Theorem 6. 


Corollary 2: Let {7,} be a second-order linear recurrence defined by 
Tn+2 r= Qlne) + aol’, - 


Then for all primes p > 3 and for all nonnegative integers b, 


Tn +2p> = QT, 4 pe + aT, (mod 2p), (17) 


where 7 is any nonnegative integer. 


Proof: Let p > 3 be a prime. By Theorem 4, congruence (17) holds modulo p for 
all n. We will show that (17) also holds modulo 2 for all ». The corollary 
will then follow since (Zz, p) = 1. 

First, suppose that Shay. Considering the characteristic polynomial f(x) 
of {7,,} modulo 2, we have 


f(#) = 2* - aye - a, = 262 - a,) (mod 2). 


Hence, the characteristic roots of {7,,} modulo 2 are PY, = ay (mod 2) and Yr, = 0 
(mod 2). As in the proof of Theorem 5, we have that if h is any nonnegative 
integer, then 
(h) (A) m 
n+ 2h Fi ay swei - ao I, (18) 


where “ and as” are defined as in equation (11). Constructing the primary 


linear recurrences {Vi} and {V6} as in the proof of Theorem 5, we observe 
that 


VO) = a, (mod 2) (19) 
for all m 2 1 and 

v2) = a, = 0 (mod 2) (20) 
for ally 21. By (12) and (18)-(20), we see that for jg = 1 or 2, 

re = yi) = a; (mod 2) (21) 


for all positive integers h. Letting h = p>, equation (18) and congruence (21) 
lead to the congruence 


Tneop» = AT 4p: + 427, (mod 2) 


which is what we wanted to show. 


Now, suppose that 2fao. Constructing the primary recurrences {yo} and 
{v2} as in the proof of Theorem 5, we see that {V{!)} and {v§2)} are each purely 
periodic modulo 2 by Lemma 2. Further, one can easily determine that the period 
of the second-order recurrence {V{!)} modulo 2 is either 2 or 3, and the period 
of the first-order recurrence {/‘4)} modulo 2 is 1. It thus follows that if we 
determine the modulus g, as in the proof of Theorem 5, then g = 2 or 3. By 
Theorem 5, if g = 2 and p is a prime such that p = 1 (mod 2), then congruence 
(17) holds modulo 2. By the corollary to Theorem 5, if g = 3 and p is a prime 
such that p = 1 or 2 (mod 3), then the congruence (17) again holds modulo 2. 
Since for any p > 3, p = 1 (mod 2) and p = 1 or 2 (mod 3), the result now fol- 
lows. 





Remark: Note that Corollary 2 to Theorem 6 generalizes Theorem 2. 
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(Submitted January 1987) 


1. Introduction 


For a real sequence {a;},.9, @9 # 0, and » € N, consider the real array 
A(k, ns u)> (kK; n) € Zx Z, which satisfies the recurrence 


A(ks ns uw) = ~ a,A(k - j,n- 1, u) (1.1) 
j20 
subject to the diagonal condition 
A(ks uks uw) = O for k > 0, (1.2) 
and the conditions 


A(O, 0, np) = 1 and A(k, n, vw) = O for k < 0. (1.3) 


We wish to use lattice path combinatorics to obtain known formulas for A(k, 
nN, u). Collectively these constitute a Lagrange inversion formula. Others have 
made similar studies; our explanations are influenced by those of Raney [18] 
and Gessel [9]. We examine specific examples of recurrences and their 
solutions, the generalized Catalan arrays. We illustrate our approach by 
enumerating certain plane trees. 


For the given sequence {a@;},.9, let a(x; u) denote 


DA, uk +1, wek, 
k20 


which we view as a diagonal sertes. In particular, let 


a(x) = a(r3; 0) = )) a,x* (the initial series), 
k20 


and let 
a(x) = a(x; 1) (the principal diagonal series). 

For any power series, let [xk] xf, 24 denote the coefficient f,. Let 
a, = [x*]a(x) = A(k, k +1, 1) [the principal diagonal of A(k, n, 1)]. 


It is immediate from (1.1) and elementary properties of formal power series 
(see [Sis f12}) that 


A(k, Ns 0) = [xkla"(x). (1.4) 
The following record solutions to (1.1, 1.2, and 1.3). 


Propositions: For m, n € Z and k, id € NW: 
n —- uk 





t.  AGés. 7s.) ARs ye Oe 2 =O. (1.5) 


2. A(ks ns u) 


LX (1-usdajyA(k - JG, n- 1, 0)3 A(k, 1, vw) = (1 - uk)a,. (1.6) 
j20 
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3. Atk, m+n+uk, rX) = DAG, mtu, MAK -Jj,nt+ulk- 7), v)- (2.7) 


Jj 20 
4. A(k, n+ uk, uw) = [xkJa"(x3 y). (1.8) 
5. (x) = a(x; p) is a unique series satisfying $(x) = a(xo"'(x)). (129) 


These are proven in Sections 3 and 4. In proving (1.5), we interpret the 
factor (nm - uk)/n. For (1.6) we interpret {A(k, 1, Wheeo- Proposition (1.7) 
is a Vandermonde-type convolution; (1.8) shows that A(k, nv, uw) is a convolution 
array; (1.9) gives a functional relationship between a(x; wu) and a(x) which 
immediately yields xa(x) as the compositional inverse of x/a(x). Correspond- 
ingly, (1.9) with (1.4) and (1.5) yields a Lagrange inversion formula; another 
is given in Section 5. 

A lattice path is a directed path in the Cartesian plane with vertices the 
lattice points (integer pairs) and with steps (directed edges) of the form ((xz, 
y), (x +u, y +v)). There will be various restrictions on (u, v); the set of 
permitted (u, v)'s is called the step set. A lattice path from (0, 0) to (k, n) 
which lies strictly above the line y = px for 0 < x < k is called a (k, n, u)- 
path. If we restrict to steps of the form ((x, y), (x + gj, y + 1)) with weight 
aj;, and if the weight of a path is the product of the weights of its steps, then 


we shall show that A(k, n, u) is the sum of the weights of the (k, n, )-paths 
for n > pk. 





FIGURE 1 
A(3, 4, 1) counts the (3, 4, 1)-paths with step set {(j, 1): 7 € N} 
anda; = 1 for j € N. [A(3, 4, 1) = C(3, 4, 1) = y3 of Example 2B).] 


2. Examples of Recurrences and Their Solutions, the Catalan Arrays 


The recurrences are defined by their initial series. A(k, n, 0) anda = 
Atk, kK +1, 1) (often in [25]) are found from (1.4) and (1.5). For reference 
Dig); e(x)s ete, BUKs ns ws Cs ms. wi, tes; and 85. 7¥,-ete. dénote the spe= 
cific a(x); A(k, ns 1), and a. Here “PA,” “CA,” and “CN” abbreviate Pascal's 
array, Catalan's array (see [21], [24]) and the Catalan numbers [11]: 1, 1, 2, 
5, 14, 42, ... . These examples are unnecessary for Sections 3 and 4. 


Example 2A: b(x) = 1+. B(k, n, 0) = () , PA, where, for € Z and ke N, 


(;.) = (n)(n - 1) «++ (n= k + 1)/k! if k > 0 and (5) x 


BK 2s A) e i se another PA, and B, * i for k =O: 
ae nai), CA (eee Table 1). 
[xk]@(x3 2)=B(k, 2k+1, 2)= weit _ CN (marked « in Table 1). 


B(w; 2) = 1 + wB2(e3 2) by (1.9). The step set {(0, 1), (1, 1)} yields (%) as 
34 
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the number of (kK, 1, 0)-paths. 


So 


n\k l 2 3 4 n\k 

















0 1 2 3 4 
5 1 a 2< -2 -3 5 1 4 9 14 14< 
4 1 2 0 -2 -1 4 1 3 5 5< 0 
3 1 l<« -1l -1 0) 3 1 Zz 2+ 0 -5 
Za 1 0 -l 0 0 2 1 1+ 0 -2 -5 
1 l<« -1l 0 0 0 1 l< 0 -1 -2 -3 
0 1 -2 2 -2 2 0 1 -1 -1 -1 -1 
-l 1 -3 5 -7 9 -1 1 -2 0 0 0 
-2 1 -4 9 -16 25 -2 l -3 2 0 0 
TABLE 1 TABLE 2 
A section of B(k, n, 2) A section of C(k, n, 1) 
Example 2B: ec(x) = >, xk = (1 - x)7!, 
k20 
. _(n+k- 1 
C(k, n, 0) = ( . ee 
C(k, n, 1), CA, are the ballot numbers [3], [16], see Table 2. 
a ‘ -(2°)>eN (marked + in Table 2). 


xy2(x) - y(x) +1=0 by (1.9). C(k, n, 9) counts the (kK, n, 0)-paths with step 
set {(0, 1), (1, 0)} with (0, 1) as the initial step. C(k, n, 0) also counts 
the (k, n, 0)-paths with step set {(j, 1): 7g € N}; see Figure l. 


Example 2C: b(x) = 1+ x’ = b(x’), where v EN, v > O. 
if if k 


B(k, n, 1) = 2 EK) <,) 


B(k, n, 0) 


vK and = 0 otherwise, a variant PA. 





Bhs 7s pO kE vK and 


QO otherwise. 


Example 2D: (x) = 1 + x2. 





C(k, n, 0) = (x72) for k even and = 0 otherwise, a variant PA. 
S apeb eg Ryd ; 
lie "Ee ( k/2 ) = 1, 0, 1, 0, 2, ..., zero-interspersed CN. 


Example 2E: G(x) = 1+ 2x + x2, 


C(k, n, 0) = () a PA with every other row missing. 


C(k, ns 1) = B(k, 2n, 2). 


= 1 2k 2 : : 
= £ + en i i 36 
Yx Bk, 2k Bs 2) k -( k h CN with first entry missing. 


Note that [x*]7¥(a) = [x*]@2(x3 2) = [x*]a7!(B(x;3 2) - 1). 


Example 2F: mx) =1+2 + x. 
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‘ , : ae! _yg(nyn titel 
M(k, n, 0) = (21 - #8)(1 - 2)7} eee wan 
Also ‘ 
M(k, n, 0) = [xk] ((1 + 2?) + x)” = 2, (72). 


t,g:n-j+2i=k 


i, = Mk, K+1, 1) = 1, 1, 2, 4, 9, 21, ..., named for Motzkin [17], who found 
them to count the ways of placing nonintersecting cords between k points on a 
circle. Note 











— k 2 a ee 1, go 2y4 
fig = leh (C + a?) + ayetl = po 7%, ( ; )teF-1] (1 + x2)d 
k \s 
“Zl; ~ 1)¥e-1 
and 
Ye = age tet et) tamktt act yD (4) ped tc + + 229 
kk k+1 k+1f59\ § 
k \- 
= - Us_y.6 
ZAg - sia 


See Example 7A and [4], [5], and [14]. 


3. Lattice Path Analysis for Propositions 1 and 2 


We use weighted paths with steps of the form ((x, y), (x + Jj, y + 1)), de- 
noted by <j» and assigned the weight aj, 7 2 0. ‘j1:n>% denotes an arbitrary 
path <j)><jo> +--+ <j, and Mini a;, denotes its weight. P(k, n, uy) denotes the 
set of all (k, 7“, u)-paths and |P(k, Ris u) | denotes the sum of the weights of 
the paths in P(k, n, uy). When appropriate, |A| denotes the cardinality of A. 

Since all (k, n, u)-paths pass through {(k - 7, n - 1):0 < j < k} exactly 


once, 
n 
zs I a;. 


<Jizn? © P(ksn, u) t=1 


|P(k, n, u)| 


n-1 
Il 2;,2; 
Osgsk (fyrn-p?E€P(K-Gsn-1,u) t=1 °° 


LD a;|P(k - j, n- 1, uv). 
O<jsk 


|P(O, 0; u) | = 1 and |P(k, uk, u) | = 0. Hence, |P(k, aie u) | satisfies (1.1), 
(22); and (1.3) fern 2 wk. ‘Thus, 


Proposition 6: |P(k, n, u)| = A(n, k) for n 2 yk. (3.1) 


We next determine |P(k, n, u) | by a "radiation" scheme, which extends the 
method used by Dvoretzky and Motzkin [7] on Barbier's ballot problem of 
counting the (k, n, w»)-paths with a two element step set. Grossman [13], [16] 
reformulated their technique as "penetrating analysis." See also [18]. 

Each path <j,:,> € P(k, n, 0) determines a sequence of cyclic permutations, 
each being a path in P(k, n, 0): 


SFisn?» SJaz1? = SJ2?6G3? 29+ SFn I> > 63:2? = SI3><Fy? oe + 61? 6I2?s 
cows SSgpnng? = ee +3 Get, (3.2) 
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Let p be the period of <j1:,>- [<j1:,>2] denotes the cyclic permutation 
class {(j1.,>> (J2:1>» +++> ‘dp:p-1>}, the set of distinct paths in (3.2). 
Let <j1:,> be a fixed [fixed until (3.5)] path in P(k, n, 0) of period p. 


‘FJizn? = 691><J2> +++ “jp? € Plkp/n, p, 0); 


is the initial subpath of ‘j,:,>. Each path in [<j,:,>] is the concatenation of 
n/p copies of a cyclic permutation of <j,:p». Distinguish the steps in ‘(j): p> 
by their index. Thus, each step in ‘j,;,> initiates a unique path in [<j):,)]. 
A step <j> is called a zero step if j = 0; otherwise it is called positive. Let 
J, = {%: 7; > 0 andi < p}, the index set of the positive distinguished steps of 
<j1:p?+ See Figure 2. 


<%)> > 


(0) -» 





FIGURE 2 


If wp = 1 and ¢J71:7> = <0)><09>23><0,><45><06><07> (the subscripts 
distinguish the steps), then 4={3, 5}. 23 = {2}, Z5 = {l, 4, 7}, 
Ze = {6}. <J1:7> and <jg:5> are shown with 3 subswaths of rays in- 
dicated on <jg;5>. 


In the following m 2 uk. Fix 7 € J, and consider the geometrical configura- 
tion of <j;41.;> © [<71;p?] where each step is a line segment (7 + 1 is replaced 
by 1 if « = p). Fromthe points on the last step of ‘j;4,,;>, namely <j;>, draw 
rays in the direction of the ray from (0, 0) through (-l, -y). See Figure 2. 
Since the terminal vertex of ‘j;41:;> is above or on the line y = ux, all rays 
must strike and be absorbed on the right side of <j;4,,.;» by the zero steps, 
positive steps being too inclined to be hit. By examining a triangle with 
vertices (0, 0), (0, ug - 1), and (j, uj), we see that the vertical width of 
the swath of rays from <j;> is uj; - 1. This swath can be _ partitioned 
into uj; - 1 equal parallel subswaths. Since each subswath passes between 
vertically adjacent lattice points, each subswath must irradiate the entire 
interior of a zero step. If Z; denotes the set of indices with respect to <j); p? 
of the zero steps which are irradiated by the rays from <j;>, then 


|2;| = ud; - 1. (3.3) 


We claim that the Z;, 7 € J,, are disjoint from one another. Suppose there 
is a zero step that is irradiated by both <j;> and ‘j;,>, where the zero step 
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appears earlier in, say, ‘j;-4,:7°>. But the configuration of <j;'4): ;/> shows 
that the step <j;> will shield this zero step from the irradiation of <j;->. 

Let Z. be the index set of clean (nonirradiated) zero steps in ‘j);p>. By 
considering the vertical and the horizontal dimensions of <Ji:p?> 


p= > [2] + losl + [Z| 


and nye 
ukp/n =U > ™ j= +. (uJ; ~ aE + le, | zi |Z, | € lJ |. 
and thus an _— — 
- uk 
\Z.| = p(n - uk) (3.4) 


As noted, each path in [‘j,:,»] is uniquely determined by its initial step 
which is a distinguished step of <j,:p». A path beginning with a positive step 
touches or is below y = ux by the first step. A path beginning with a zero 
step touches y = yx for the first time on its radiating positive step. Thus, 
the paths beginning with a clean zero step are precisely those belonging to 
PUR, 25 i)... “By (324), we ‘have 


Lemma: The number of paths in [<j,:,>] 9 P(k, n, yu) is 


_ p(n = pk) _ n — wk 4 
IZ. alee mumaiaings wee | (3.5) 
Since every path in a cyclic permutation class has the same weight and 
since the classes are disjoint with union P(k, n, 0), 


n-wuk,,, B 
|P(k, n, u)| = > Ea LTD a, (sum over all 
ure c.p. classes) 


A(k, ns wu) 


Po TG : n n- uk 
ge NEST Beg, = ons ts OF 


n — uk 


ul 


Atk, ny BW. 


Thus, a formula for A(k, n, ut) has been constructed for nm = uk. Simple arith- 
metic shows that this formula satisfies (1.1), (1.2), and (1.3) for n # 0 and 
k 2 03 hence (1.5) is valid. 

A second realization of the contribution of each cyclic permutation class 
to | P(k, Ny u)| establishes (1.6). By (3.3) and (3.5), the weight contributed 
by [<j1:,2] for n 2 pk is 


iol fas, = OO - wieday, May, +E 4] EG - weds, May] 


tEZe GEd, sé Z, 
(since jz, = jg = 0 for zero steps and 


+ (1 - U4 )az, Nes. since the term in { } is 0) 


Pp 
a i= Wd taj, IT a;,. 
t=1 tet 


Summing over all cyclic permutation classes yields 


(sum over 


n 
A(k, n, u) = |P(k, n, u)| = p> (1 - wd )a;, fas. all paths) 


Fiz EP(k, n, 0) 
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n 
y+ ( - ud a; : aj, [Since all paths pass through the line 
j 20 ‘Jain? t= 2 y = 1, the second sum is over all paths 
from ¢7, 1) t6 (kh, #).] 


> (1 - uja;|P(k ~ 35? = Ls 0) | (weights are transition invariant) 
j20 


> (1 - uj )ajzA(k = ds n-1, 0). 
j20 


Thus we have constructed the formula of (1.6) for m = uk. If A(k, mn) momen- 
tarily denotes this formula, then it is easily shown that A(k, n) satisfies (1.1) 
for (k,; n) € Nx Z. Since A(k, n) and A(k, n, u) agree when nv > uk, they agree 
for all (k, n) € N x Z, yielding (1.6). 

Equation (1.6) yields a nice interpretation for A(k, n, u) on both sides of 
y= yx andn 21. Retaining the definitions of this section, reassign the 
weight of (1 - uj)aj to the initial steps ((0, 0), (j, 1)), j 2 0. Then A(k, 
n, uw) is the sum of the modified weights of all unrestricted paths from (0, 0) 
to Gi, Kk). 

With a’(x) denoting the usual formal derivative, immediately (1.6) is equi- 
valent to (similar to a result in [1]) 


A(k, n, vw) = [x*] (a(x) - uxa'(x))a"-1(x%) for ne Z. (3.6) 


4. The Proofs of Propositions 3, 4, and 5 


We establish (1.7), a useful generalized Vandermonde-type convolution [10], 
[16]. Then, using the tractable notation of series, we reformulate both the 
convolution and (1.1) in terms of diagonal series. 

First we give a lattice path proof of (1.7) for m, n 2 0 and m + uk 2 dk. 
Since ((x, y), (x+j, y+1)), jg 2 0, is the form of the lattice steps, any path 
in P(k, m+n+ uk, X) must intersect the line M = {(j,m+t+uj):0< j < k}. 
Since the weight of a path is invariant under translation, the sum of the 
weights of the paths from (j, m+ uj) to (kK, m+n + uk) which remain above Y/Y is 
|P(k-j,n+u(k- 7), u)|. Hence the sum of the weights of the paths in P(k, 
m+n+uk, X) that pass through M for a last time at (Jj, m+ uj) is the product 


IP(j, m+ ug, HD||PRK-gG,n+uk- 7), w)]- 


Summing over M and putting A(x, y, u) = |P(x, Y> uv) | yields (1.7) in this case. 
Now for meé Z and n > 0,'(1.7) can be proved by induction on the value of 
n+ uk by observing that, for n + pt - 1 <n + uk, 


> aA(kK -t,m+n+ ut -1+u(k -7), A) 


120 

= YaAj, mtug, MAK -i-j,ntui-l+ulk-7t- 9), w) 
~20 j20 

= DA(j, m+ ug, A) DX aACk-groi,n-1l+uk- 7), u)- 
J20 ~20 


The case for me Z and n < 0 can be proved by induction with respect to -m” upon 
noting that (1.1) yields 


A(k, me Or = by + uk, LW) 


=A(k, m+n+uk, u) - YajsA(kK-G,m+ n-1+ug) + ulk - J)> wd) 
gj 20 


and that -2 + 1- uj < -n. 
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Equivalent to (1.8) is 


a(x; un) = DL A(k, uk +n, w)xk for ne Z. 
k20 


For n 2 0, this can be proved inductively since, by (1.7), 


a"*l(es yp) = ales wa"(e3 vu) = DACk, uk + 1, uak DAK, uk +n, wack 
k20 k20 


= DA, pukK +n4+1, v)axk. 
k20 


The case for n 2 0 and (1.7) yields 


ange) 


a-"(a) DACk, uk +n, u)ak DO Ak, uk - n, wack 
k20 k20 


a-"(x)a"(a) PACK, uk - nm, weak = YAK, uk - n, u)ak. 
k20 k20 


As in [9], equation (1.8) has the following meaning for nm = wk: Since each 
(k, uk +n, w)-path must sequentially intersect for a last time each of the lines 
y=ucx+d for 1<t<n, each(k, uk +n, )-path is an n-fold concatenation of 
(k, ug + 1, w)-paths for various j. Correspondingly, the total weight of the 
(k, uk +n, w)-paths is a coefficient of an n-fold convolution of a(x; u). 

Moreover, since each (k, wk + 1, y)-path intersects the line y =wpk only 
preceding its last step, the set of (k, uk + 1, u)-paths is the disjoint union 


k 
P(k, uk +1, uw) = U {51:4 67? t Fis pe? is a 
as ‘ : ‘ 
k (k i ds u(k Zz a + Ud >» u)-path} 
=U {<91:44F? ? SJ: u,? is a uj-fold concatenation of 
j=0 


d ; , 
various (j, uj + 1, u)-paths}. 


More precisely, we have that 


a(x; u) = ) A(k, uk +1, udack 
k20 


=D Lawak - 9, wk - f) + ud, waked 
k20 720 


= Da,edVak - fs ulk - §) + uj, wak-d 
j20 k20 


= Yajxia" (x3 u) = 2 az(wo4(x3 u))’. 
j20 j20 


This establishes (1.9) since comparing coefficients shows the uniqueness. 
As a consequence of (1.9), we have 


Proposition 7: For each pe WN, if a(@; yp) is taken as the initial series, then 
a(x; u + 1) is the corresponding principal diagonal series. 
Proof: If a(x) denotes the principal diagonal for a(x; u), then by (1.9), 

G(x) = a(wa(x); u) = a(xa(x)[a(ea(x); u)]") 


a(xa(x) (a(x) ]4) = a(x[a(ae)]#*!). 


But a(x) must be a(x; u + 1) by the uniqueness in (1.9). 
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5. A Lagrange Inversion Formula 
A common Lagrange inversion formula [3], [12], is included as it easily 


follows (1.4), (1.5), (1.6), and (1.9). See [3], [8}, (97; [12]. and [18) for 
more general formulas. 


Proposition 8: For any initial series a(x), there exists a unique series 


w(x) = }) w,xk 
k21 


such that w(x) =xa(w(x)). Moreover, if f(x) is a formal Laurent series so that 


f(x) => f,«* for some t € Z, 
k2t 


[x"]f(w(a)) = qlan-1]f"(x)a"(z) for n = 0, 
[x9] f(a) + [x7] f'(x)log(a(x)a~!(0)) for n = 0. 


Proof: By (1.9), w(x) = xa(x) is the unique solution. It suffices to show the 
second part for f(x) = xk, ke Z. Forn= 0, 


[2"] (wa(ae))* = [2"-k]qk (x) = An - ky ny 1) = (a"-E at(a) 
a dat} ]kek-1e"(a). 


As noted in [12], 


0 a1 }L (wk10g (a(x)a71(0))) 
dx 


= [a7!]kak-llog(a(x)a7!(0)) + [x7 Jeka'(x)a7} (x). 


For n = 0, it follows from (3.6) that 
[x ]Jw*(x) = [x2 Jakak(x) = [x7 ]ok(a) = A(-k, 0, 1) 
= [x*](1 - xa'(x)a7“!(x)) = [x7 ]kxk-11log(a(x)a-!(0)). 


6. More Examples of Recurrences 


Example 6A: r(x) = 1+ (wt+1) Do xk = (1 + wr) (1 - x)7!. 
k21 


R(k, n, 0) = [x*] >> (% ode yo ‘ or ant 


j20 i20 . 


ee. 


i20 ~" t 
Note how Yr, relates to a, of Section 3. Also, 


R(k, Ns u) = R(k, 1 = Bs uw) + wR(k - Ls mets; u) + R(k - 1, Ms w)> 


and R(k, n, uw) is the sum of the weights of the (k, n, u)-paths with the step 
set {(0, 1), (1, 1), (1, 0)} where (1,1) has weight w. As shown in [21], [22], 
or from R(k, n, 0), we have 


Rik, ns 1) = D wi(" * an Vock «Js ga Ds 
j20 J 
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For w = 1, see Table 3, where 


. x= J 
P, %, | a Mas 
the r-Schréder numbers, are marked. See [19], [21], and [22]. Note that 


p, = 20, (k > 0)3 


see Example 6E. xp*(x) + (we - 1)p(x) + 1 = 0 by (1.9). 


n\k 0 1 2 3 4 n\k 0 1 Pe 3 4 
4 1 6 16 22+ 0 4 1 3 7 ll< 0 
_ 1 4 6< 0 -22 3 1 Z 3< 0 -22 
Zz 1 2+ 0 -6 -16 2 1 l< 0 -6 -28 
1 1< 0 -2 -4 -6 1 l< 0 -2 -8 -24 
0 1 -2 0 -2 0 0 1 -1 -3 -7 -15 

-1 bt -4 6 -8 10 -l 1 -2 -3 -4 -5 

TABLE 3 TABLE 4 
RER, ns 1) foer-w = 2 SCk; 25. 1). for w= 1 


General Example 6B: Given any sequence wW ;, W9, W3, ...-, consider the step set 
{(0, 1)} U{(j, 0): 7 > 0} where (0, 1) has weight 1 and (j, 0) has weight w;. 
If A(k, n, u) is the sum of the weights of the (k, n, y)-paths (the initial step 
must be vertical), we have 


A(k, nN, U) = A(k, n= ts. st D w,Atk = Js Ns Hs 
j>0 


It follows inductively that this A(k, n, u) satisfies (1.1) for f{a,},,9 defined 
by ay = 1, @, = Wis Ag = Wo + WW): Az = Wz + WoW) + WW. + WWW], and in gen- 
eral 


—* me .. I] ¥;,- 
typtigt---+t1,2#k letsm 
Hence, 


ae) < Fae = 14+ | ( Zz, w,0t) = (2 -> v,et)-), See [26]. 
k2=0 k2z1\te2l t21 
Example 6C: In Example 6B put w,] = W2= 1 and w; = 0 for 7 > 2. 
a(x) = (1-2 - x2)7!, 


Thus, a, = 1, 1, 2, 3, 5, ..., the Fibonacci numbers. See 7B. 


Example 6D: In Example 6B put w; = 1 for 7 = v and = O otherwise. 
G(x) = (1 + 2’)! = e(e’) of (2.2). 


Pia 1 
¥, 20k, k+1, 1) = CK k +1, v) = Ga4( 





ee a 


for k = vK and = 0 otherwise. 
Example 6E: In Example 6B put w; = w for 7% > 0: 


s(x) = 1+ Y Ww t+ 1)k-1ek = (1 - we) (1 - (vw + 1)z)7!. 
k>0 
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S(k, n, 0) =F cpivie owe ec ; i 


When w = 1, see Table 4 for S(k, n, 1) where 0, =S(k, kK + 1, 1), the s-Schréder 
numbers [23], are marked; see Example 7B. 

From (w + 1)xo2(x) - (1 + wr)o(x) +1 = 0 [by (1.9)] and the last identity 
of Example 6A, one can show: 

(i) (wt 1)(o(x) - 1) = p(x) - 1, 
(ii) (1 + weo(x)) = (1 - weo(x))-!, and 

(iii) o(x) = (1 - xp(x))7!. 

These are illustrated in Tables 3 and 4: (i) relates the principal diagonals. 


(ii) and (iii) relate the partial row sums in the triangle above the zeros in 


one array to the principal diagonal in the other, as generalized in the follow- 
ing: 


n-1 
Proposition 9: If t,= }. w"-*A(k, n, 1) = ¥ wiA(n- fj, n, 1), a weighted par- 
k=0 j 20 


tial row sum, then t(x) = )> t,x2" = (1 - wea(x))7!. 
n=0 


Proof: t(x) = > > wiA(n - jsn, lar= > wid > An-g,n-g+ds 1)"77 
n20 j20 j20 n2=0 

Y wiedad (x). 

J20 


This extends a result in [20]. 


7. Enumerating Plane Trees 


Informally, a rooted plane tree is an unlabeled tree which is oriented in 
the plane so that it branches upward from a root (a distinguished vertex which 
need not be univalent) to the leaves. Two plane trees are equal if one can be 
continuously transformed into the other in the plane so that the nonroot 
vertices remain above the level of the root. A more formal definition is given 


by Klarner [14] and [16]. A planted plane tree is a plane tree with univalent 
root. See Figure 3. 


yy Y y 


FIGURE 3 
This illustrates the planted plane trees with 4 edges and no 
degree restriction. These trees are listed as they correspond 
to the paths of Figure 1 under Bijection A. MThe numbers 
indicate the order of growth. 


Here we enumerate rooted plane trees with vertex degree restrictions by 
establishing bijections between the trees and previously counted lattice paths. 
Equivalently, one can establish directly a recurrence for the tree counts in 
the form of (1.1). The following examples are enumerated by other methods in 
[4], [5], [7], [12], [14], [15], [16], and [22]. One common method for planted 
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plane trees is to establish a functional equation for a generating function and 


then solve the equation perhaps, and not surprisingly, by the Lagrange 
inversion formula as in [12]. 


Bijection A: Counting rooted plane trees with respect to the number of edges 


Let T(n, p) be the set of such trees with n edges and root degree p. Let D 
be a specified set of permitted degrees for the nonroot vertices. Let P(n - 0, 
n3;D) be the set of (n - 9, n, 1)-paths with step set 


{(j, 1): j + 1 € D}. 





A bijection from P(nm - p, n; D) to T(n, p) is defined inductively. The trivial 
zero-length path in P(0, 0; D) corresponds to the tree consisting of just a 
root. A path in P(n-p, 3D) and the corresponding tree in 7(n, p) can be 
extended in two ways. (i) The path can be extended to a path in P(n-p,"+ 
1; D) by attaching a new step (0, 1), while the corresponding tree is extended 
to a tree in T(n +1, p+ 1) by grafting a new left-most edge to the root. 
Such a new step corresponds to a new leaf. (ii) For 7 +1€D and j <p, the 
path can be extended to a path in P(n -p + J, n + 1; D) by attaching the step 
(j, 1) while the corresponding tree is extended to a tree in Tin+1,p- 7 + 1) 
by cutting at the root the j left-most incident edges and then grafting the 
lower vertices of these edges to the upper vertex of a new left-most edge 
incident to the root. Thus, a (j, 1) step corresponds to a new vertex of 
degree j + 1. Hence, 
IZ7@, )| =A -0, n, 1) = Stare Fad)", 


j+1eD 





Bijection B: Counting rooted plane trees with respect to the number of leaves 


Modify the scheme of Bijection A by replacing the underlined phrases sequen- 
tially by: 7 leaves; let D (2:¢ D); 10; D3 vutgG. O):7 + 2eD!: g+2eED 
andj <p - 1; Pin-p +f, n3 D)3; step (Jj, 0); Tim, p - j)3 the 7 +13 (J, 0) 
step; degree 7 + 2. Thus, by Example 6B with w; =.1 if 7+ 2€Dand = 0. other- 
wise, seeien 

IT, p)| = A(n - p, nm, 1) = BiatPH( ~ + 4) ° 
g+2ED 


Example 7A: Applications of Bijection A 


If D=N- {0}, no degree restriction, |7T(n, 1)| = C(n - 1, n, 1) = Yn-13 
see Example 2B and Figure 3. For D = {1, 3}, trivalent planted trees, |T(n, 1) | 
=C(n-1, n, 1) = ¥,-1; of 2D. For D = {1, v + 1}, use 2C. For D= {1, 2, 3}, 
no vertex has degree greater than 3, |T(n, 1) | =M,)(n - 1, n, 1)=Uy-13 see 2F. 

If D=N - {0, 2}, no bivalent nonroot vertices, let d(x) = (1 - x)7! - «x. 


[T(n, p)| = plan ]d"(x) "i n rn i : p ): 
iri, | =D Cor*(F Teo: 


(x + 1)6(x) = (1 - x6(x))7} by (1.9). One can show 
(2 + 1)(6(x) - 1) = x(i(x) - 1) (thus, 6,-)+ 6, = Uy-1], 2 > 0). 
Therefore, by Proposition 9, 


n-1 
ee Dtks ms 1) = tin-1- 
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Hence 
Zz |T(, 0) | = eas 
p21 

as in [4]. 


Example 7B: Applications of Bijection B 


An immediate source for such trees is the problem [3] of counting the ways 
to bracket m nonassociative, noncommutative factors so that the number of fac- 
tors associated by a pair of brackets is restricted to some set B. If B = {2}, 
we have the problem of Catalan [2], 1838. There is a simple bijection between 


the usual pairwise bracketings on m factors and the planted plane trees with D 
= {1, 3} and m leaves. For D = {1, 3}, 


|T(n, p)| = ore Pher - x)". 


|T(n, 1) | = Yn-1,» the appropriately named sequence of Example 2B. 


If B =N - {0, 1}, we have the problem of Schréder [23], 1870. If n= 4, 
the bracketings are 


(a(b(ed)), (a((be)d)), (a(bed)), (((ab)e)d), ((a(be))d), 
((abe)d), ((ab)(ed)), (a(be)d), (ab(ed)), ((ab)ed), (abed). 


There is a simple bijection between the unrestricted bracketings on n factors 
and the planted trees with m leaves and no bivalent vertices. For D=N - {0, 
2}, refer to Example 6E with w = 1: 


|T(n, 0)| = S(n- 9, 7, 1) and |7(m, 1)| = o,-}. 
If D = {1, v + 2}, 


|7T(m, 1) |= C(m - 1, 2, 1) = Yn-13 
refer to Example 6D. If D = {1, 3, 4}, refer to Example 6C. 
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It is well known that the powers of the matrix 


_ ? 
[ll 2] 


are matrices of the form 
. | an i 
Fue Pn+2 
whose entries are the Fibonacci numbers. If these matrices are normalized by 


dividing by ? , so that the first entry is always 1, then the resulting se- 
quence converges to 


. 


f2 
where f = (1 + ¥5)/2 is the golden ratio. 
Moore [2] noticed that if the 2 in M is replaced by 1 + x and the same pro- 
cedure (taking powers of M and normalizing to obtain a 1 in the first entry) is 
performed then the resulting sequence seems to converge to a matrix of the same 


form with f = (x + vx2 +4)/2. These observations naturally suggest the fol- 
lowing questions. 


(1) 


1. If we start with any symmetric 2 x 2 matrix M, with positive integral 


entries, does a similar phenomenon occur and, if so, what is the corresponding 
value of f? 


2. If there is convergence with f = (a + Vd) /b » what are the values of d 
that can occur, i.e., in what quadratic number fields do we find such f? 


Since we normalize at each step, we can assume that 
t y 
M= P 
l+ez 


Diagonalizing M gives M = PDP~! with 


1 1 (x + 2 + vd)/2 0 
P= D= 


> 


(2 + vd)/2y (x - Vd)/2y 0) (x + 2 - vd)/2 


where d = x2 + 4y2. Thus we have M" = PD"P-!, When we normalize D to make the 


leading entry 1, the second diagonal entry is less than one and so the sequence 
of its powers converges to 


E 6G 


0 60 
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and so PD P~!, when normalized, can be seen to converge to a matrix of the form 
(1) with f = (x + Vd)/2y, where d = x* + 4y2. 

The set cf d that can be written as x* + 4y2 can be easily found, since we 
know what numbers can be written as the sum of two squares. (A positive 
integer is the sum of two squares if, when factored, all its prime factors con- 
gruent to 3 modulo 4 occur with even exponent, see, e.g., [1].) If d is odd, 
then d is the sum of two squares if and only if it is of the form x? + (2y) 2 = 
x? + 4y2, since one of the terms must be even. If d is even, then d = x2 + 4y2 
if and only if d = 4m, where m can be written as the sum of two squares. Since 
dis even, x is even, and sodis divisible by 4. The 4 can then be factored 
out giving m as the sum of two squares. The converse is also easy. Thus, d is 
of the form x2 + 4y2 exactly if, when factored, all its odd prime factors 
congruent to 3 modulo 4 occur with even exponent and 2 does not occur with ex- 
ponent l. 
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Introduction 


The three related classical series for representing real numbers as the 
sums of reciprocals of integers were all studied toward the end of the nine- 
teenth century. These are, respectively, the series of Sylvester, Engel, and 
Liiroth (see Perron [1]). More precisely, given any real number A there exist 
three (different) sequences of integers {a;} such that 

1 1 1 


a, +—t+—t— tee 


(i) A 
Q@ @ 44% 


where a, 2 a5 Qz4 1 2 a;(a; =— b+ 1 for -¢ > 1, 
1 1 1 


(ii) A + + 
Q@, &)4g 25a 





Nt 
& 
Oo 
= 
| 
as 


where a, > 2; a; 2 a, for 4s hs 


1 1 1 
C—— sees 


1 
iii) A=a, + — + ——— > — + ——_————_—_—_ 
(144) 0 a, (a,-1)a, a, (a, - 1)a,(a, - l)a, ay 


where a; 2 2 for i 2 1. 


Observe that as we move from the Sylvester series (i) to the Liiroth series 
(iii), the denominators in the expansion become increasingly more complex while 
at the same time the growth conditions on the digits a; become simpler. We now 
generalize the expansions in (i) and (ii) above, to obtain new representations 
for real numbers that depend on a power kK > 0. These new representations have 
the desirable property of having terms only slightly more complex than in (i) 
and (ii) above, yet their digits need satisfy only mild growth conditions. Two 
different sets of algorithms leading to results of the types mentioned are con- 
sidered. We state the main results in the case where the digits a; grow least. 


Given any fixed real k 2 1 and any real number A, there exist sequences of 
integers {a;} such that 


(i) A=a,+ ho ae et 
i = ee 
ns ee 


where a;,, 2 @; 2 2 for 7 2 1, and for 7 sufficiently large, 


t 
a; +1<s a3, s ilk, * hs 


; 1 1 1 
CHE) A 9 hy 6 ee Fe + «05 
ay (a4) (414743) 


where a) 


ul 
No 
. 
_ 
IA 


a,s 2 for ~ = 2, and a; = 2 infinitely often, 
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+ : + 2 + : + 
a es freee ae ee 
0 k k 

ay a;a, (a,a,) a3 





(iii) A 


where {- 7 A a a, < 2 fort. =. 2;end a; = 2 infinitely often. 


Since only the digits 1 and 2 are used, the representations (ii) and (iii) 
above could be regarded as being analogous in some fashion to the binary repre- 
sentation for real numbers. Expansions of the above form where the digits have 
no upper bounds are also considered. We note in particular that, by setting k 
= 1 in the above results, we obtain expansions for real numbers with the same 
form as the Sylvester and Engel series but whose digits are considerably 
smaller. In addition, when k is a positive integer, rattonal numbers have 
representations of types (ii) and (iii) above for which the digits a; become 
pertodic. This condition is analogous to that of the Luroth series when A is 
rational. 

The paper is set out as follows. In Section 2, we consider kt! power ana- 
logues of the Sylvester series. In Section 3, we consider kt power analogues 
of the Engel series. Finally, in Section 4, we consider kt power expansions 
that are related to a simplified version of the Liiroth series. 

For convenience we introduce the following notational conventions. The 
lower case letters a; and a, denote integers throughout the paper. Further- 


more, unless otherwise stated, the lower case letter k represents a positive 
real number. 


2. Generalizations of Sylvester Series 
We introduce two different algorithms that lead to a kt" power generaliza- 
tion of the series of Sylvester. The first coincides with the ordinary Sylves- 


ter algorithm for kK = 1. The second leads to a restricted growth of the digits 
in all cases, including k = l. 


Theorem 2.1: Let k > 0. Every real number A has a representation 


A aed Cones 
=@a — —_— —_= eee 
0 k k k ’ 
a; @y 3% 


where: 


if kK > 1, thena;,, 2 a; 2 2 for t 2 1, and for 7 sufficiently large, 
Qe4] 2 a; + Es 


dt Ox? k “<1, "then: Q@.4, 2 Ogle; - 1) 4+ 1 for t 21, a, 22. 


Lt 
Proof: In order to obtain this result, we introduce the following algorithm: 
Given any real number A, let A, =A- ap» 0 < Ay Sb 


Then we recursively define 





1 
a, = ae +1 forn & 1, 4, > S, 
n 


where 


A.., =4, - 


for A, > 0. 


n+l 


Q | — 
x= 


First, repeated application of the algorithm yields 


1 1 1 1 
A= a) +A, =a) t—etAy = eee wma, +t —mt—awtere t— +A. 
0 1 0 ak 2 0 ak ak ak ntl 
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Now a, = [1/A}/*] + 1 implies that, for 0 < A,, (a, - 1)k < 1/A, < ak. Thus, 


7 n 


and provided a, = 2 (0 <A, < 1) 
A < sie tana 
" (a, - 1k 


Now 0 < A, < 1 implies a, 2 2 and 


> 


awe 1 1 1 0 
oe” oe a 
| . 5 


Continuing this process inductively we see that A, > 0 for all m. Furthermore, 
since {A,} is a strictly decreasing sequence of positive values, we deduce that 
Any] 2 Ay 2 2 for n => 1. Therefore, 





‘ 4 Ps 1 1 ak- (a, - 1) 
a Paes 2 ane a kk 
b a Gy) ax (a, - 1)*ap 
Thus, 
~ 1)kak 
ak, > Be 4 ntotanns, 
Ane] an ~ a, - 1) 


In the case 0 < k < I, ak - (a, - 1)k < 1; so 
Gigs 2 (ey ~ Tay + 1s 0 2 1. 
In the case k > 1, we have 
An+) 2 an + 1 
provided 
(ay - 1)* 
ak - (a, - 1)* 


This is true if a, 2 21/k (Qik 1). and if A, < (21/k - 1)*/2. On the contrary, 
suppose that 


gi/k 
45 a,-s [sz— |- e(k), 
say. Then 


1 1 


A aA, - < A, - —. 
=e k (e(k))* 


a 


Now, either A,,, < (2!/* - 1)*/2 or 





A =A ; < alias fa < e(k) 
— 7“ exe. "— (e(k))k ha 


Thus, at each stage, A,4; is decreasing by at least a fixed COMME MNEs OS after a 
finite number of steps we must reach a stage at which A; < (21/k - 1)*/2. The 
result for k > 1 now follows, since 


1 1 
As < 


<*> 0 as rn rr, 
eee Gun DO OM) +8 = 1)" 


For 0< kK < 1, a4) 2 (an - l)a, + 1 2 a, + 1 as a, 2 23 hence, 
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1 
A, < m+b* >O0O asn>o, 
n 


A slight modification to the algorithm leads to the following results. 


Theorem 2.2: Let k > 0. Every real number A has a representation 


where 
et 1 2a, = 2 for an 2 t, 


and for n sufficiently large, 


OS tg, 52", 1, 


Proof: We use the same algorithm as previously, except that now we let 








[(7)""] lf 5&3 
a, =|(— + er 2 225 > U. 
n he n 
As before, a,,, 2 a 2 2 for n 2 1, but now 
2 
—= <A, < ——, if 0<A, <1 
aU Ge <1} . 
Therefore, for n 2 l, 
1 1 
A aA eek > ee 
n+l n ak ak 
and 
2 DTK 
Gna, 21+ la) | < 1+ 2ll/kq, 
Also, 
eck 2 1 _ 2ak - (a, - 1)* 
ee ie ak(a, - 1k 
So 
2 2ak(a, - 1)* 
ae? oa ee ee 
n+l 2a,, — (4, — 1) 


and we have a,,, 2 a, + 1, provided that 


2(a, - 1)* 
2ak - (c, - 1)* 


This is easily seen to be true if a, 2 3l/k7(3l/k _ 21/k) and if 


2(31/k es gl/k yk 


ay < % 


An argument similar to that used in the previous proof shows that these 
conditions must hold after at most a finite number of steps. Thereafter, as 
before, A, > 0 as n+ ™, and the result follows. 


We note in particular that by setting kK = 1 in Theorem 2.2 we get an ana- 
logue of the Sylvester series 
1 1 1 


A= aq + x + a +. a3 + covey 


where a,,, 2 a, 2 3 for n 2 1 and for n sufficiently large 


a, + 1 < Gye, S 2a, + 1. 
52 
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for nm 2 1 of Sylvester. However, under these weaker conditions we no longer 
obtain uniqueness for the expansions. For example, if we let instead 


a, = (zy) +1 


where m >1 is a fixed constant, we obtain a new expansion for A where the 
digits satisfy very similar growth conditions. 

As a particular case of these expansions, we note that, by definition, the 
Riemann zeta function ¢(k) for k > 1 has expansion 


1 1 1 
c(k) = 1+ ok + 3k + ae 


This is a much milder growth condition than the condition a,4; = a,(a, - 1) + 1 


+ ccs 


Euler's well-known formula for [(2m), m= 1, 2, 3, ..., then yields 


2m 
2m -1 : = a. i Je 
ee Comyy  ! * Gem + Zam tam Foote 





where B, is a Bernoulli number. 


3. A Generalization of the Engel Series 


Using algorithms essentially similar to those introduced in Section 2, we 
obtain kt® power analogues of the Engel series. 


Theorem 3.1: Let k > 0. Every real number A has a representation 


1 1 1 
Be Gg yO ere hn #34 
’ ay (a,a,)* (a,a,a3)* 
where: 


ff kK > 1, then a, 2 2, a; 2 1 for t 21, and a; = 2 infinitely often, 


41f 0 < k < 1,. then @;4, 2 4; 2 2 for pk 


Proof: We make use of the following algorithm. Given any real number A, let 4, 
=A- a), 0< A, <1. Then we recursively define 





1 
a, = aan | + ) for 2 > 13.4, > 0, 
where _ 


Ans) = An - 1 for A, > 0. 


First, repeated application of the above algorithm yields 


ee 
Aza, +A, =a, +—j += °° 
a a 
1 1 
1 1 1 A 
© (0 ee ee eee n+l i 
a; (aay) (a,a «++ Ay) (aa +++ Ay) 


Again, if A, > 0, we have Ak > l/a,, and if also a, 2 2, then 


paresis 
(a, - 1) 
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Now 0 < A, < 1 implies a, 2 2 and Ao =afA - 1> 0, thus a 21. Continuing the 
process inductively we see that A, > 0 and hence a, 2 1 for all m. We consider 
now the case kK > 1. Suppose a, 2 2, then 


k 1 k 
on -le= ( + ——— -1< 2k -1, 
ay 1 


k 
A =a,A, -1 ee 
n+l nn ‘= 1)k 


lA 


since we have assumed a, = 2. Now, if A,+,; < 1, then a,4,; 2 2. Otherwise ay+) 


= : and A,42 = An+] - 1. Continuing this process, we see that after at most 
[2* - 1] steps with 
Qan+i = 1k =A 


we must reach a stage at which 


nt+itl mera ee _ 


An+j S 1 and ayy; 2 2. 


We deduce that the sequence {A,,} is bounded above by 2-1 for all m. Further- 
more, there exists a sequence of integers n, = 0 < mp < n3 < ..- such that 


0<A S152 Sei 


n; +1 n; +1 
and a, = 1 for all other nm > 1. Then 


A 
o< n; +1 = 1 


ee ee ae a 
Kinase Fre 


and so S,. + A as t > ~, where 


1 1 1 A 
Sy = ay += +—— + + + ———__ =A - _ 


k k k : 
a; (a,a,) (a; +++ ay) (a, «+. a,)* 


Now let n;-; <n <n;. Then Sy, , < Sn <Sn,, andn+e iff t+. SoS, +A 


as n> », i.e., the series converges. For the case 0 < k < 1, if a, 2 2, then 





k 
1 (a, = 1)* a 
k Pe n = n 
ne 2 ness, COE EQ. 5 eee = hi 

r Ay “45 ~<dte,> (a, - 1)* 
Now for k < 1, ak - (a, - 1)* < 1 and since a, 2 2 we deduce that a,,, 2 a, 2 2 
for:all n 21. Thus, 

Ans 1 





Speer : > OAs, > o, 
(a, «++ Gy) . 
and again the series converges. 
This result gives the ordinary Engel series for k = l. 
It is possible to further restrict the growth of the digits a;, so that for 


t 2 1 they need only take on the values a; =1 and a; = 2, for any k 2 l. 


Theorem 3.2: Let k > 0. Every real number A has a representation 


1 1 1 
A =a) +— +——_ 


k oo ee ee es 
a; (a4) (414943) 


where: 


Ef kK = 1, then a, = 2, 1 < a; < 2 for t 2 2, and a; = 2 infinitely often, 


54 


[Feb. 





REPRESENTATIONS FOR REAL NUMBERS VIA kt POWERS OF INTEGERS 
If 0<k< 1, then a, 2 2, 1s a, < 1+ [(2)/*] for ¢ 21, 
and a; 2 2 infinitely often. 


Proof: We use the same algorithm as in Theorem 3.1 except that now we let 


Ay +4 «digg 14 Ay = 2, 


2 \1/k 
a, = (2) | +f torn 2 tA, > 0. 


and 


As in the previous result, A, > 0 and Gd, = Iter ali w 2 1. Also; 


An >= 
n 


which implies 
k 
Ag, * GpAy ~ 1 > 1 
and, in the case k 2 1, 


2 ybix 
ayer 71+ | | eae 


n+] 





Thus, 1 < a, < 2 for n 2 2, and (since 1 < A, < 2) a, 2 2. Also provided a, = 
(the case a, > 2 cannot occur for n 2 1, by the preceding inequalities) we get 
2 2a 
A/S ——eeaecaa me 2 and A,,, $ ars 

(a, - 1) (a, - 1) 


since we assumed a, = 2. Now, in the same way as in the previous theorem, after 
at most [2**1] - 1] steps of a,4; = 1, we must reach a stage at which An+z <2 
and a,4; = 2. Therefore, the sequence {A,} is bounded above by 2**! - 1 for all 
n 21. The convergence of the series for kK = 1 is now shown in exactly the same 
way as in the previous theorem. The proof for 0 < k < lis exactly the same ex- 
cept that a, 2 2 and forn21, 


l <a, <1 + [21/4]. 


-1= 2ktl -], 


In particular, by setting kK = 1 in Theorem 3.2, we get an analogue of the 
Engel series 
1 1 1 


A=a,+—+ ee eee 
Q@, @)4, &)25a3 





where a; = 2, 1 < a; < 2 for ~ 2 2 and a; = 2 infinitely often. Compare this 


to the growth condition a;4; 2a; 2 2 of Engel. Again under these weaker 
conditions the expansion is not unique. For example, in Section 4 we consider 
a different algorithm which for k = 1 gives another series with the same form 


and conditions on the digits, as the series noted here. 

We note as well that, if we had defined A; = A - ag with 0 < A; < 1, as we 
did in Theorem 3.1, the digits obtained would have satisfied the same con- 
ditions as above for 7 2 2, but would have had a, > 2 if 0 < A, < 1/2*k-1, The 
representation thus obtained would no longer be entirely in a "binary" form. 

The representation of rational numbers when kK takes on integer values l, 2, 
3, ... is also of interest. The condition that holds, i.e., that A is rational 
if and only if the digits in the expansion eventually become periodic, 
corresponds to the criterion for the representation of rational numbers via the 
Liiroth series. The result below applies to both the algorithms of Theorem 3.1 
and Theorem 3.2. 
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Proposition 3.3: Let k = 1, 2, 3, ... . The digits in the kt" power expansions 
of Theorem 3.1 (or Theorem 3.2) become periodic if and only if A is rational. 


Proof: Suppose firstly that A, = p/q is rational (with p, q «€ N). Then, since 
k € N, each A, is also rational, with 


hg ok + 1 eek a A = 4 
oe aka, +b =, 
where a€ N, DE Z. Now, for the first algorithm (Theorem 3.1) we have 
Pp 
0<A, = - < 2k -1. 
Thus, every 
sz (2k - 1)q 
Ane wae ae ae vig oe 
qq q q 


and so there exist m, n € N such that A, = A,4,. Then the algorithm applied to 
An+m Zives the same successive digits as when applied to A,, i.e., the digits 
become periodic. The same argument applies in the case of the second algorithm 
except that now 0 < A, = p,/q < 2*t+! - 1. 

Conversely, suppose that eventually a, = Qy4m- If we use the notation 


X pe ; + + : 
ee ee oe ee ee ee ee 
ay = (aap) (2). sss Mya? 
and let a, = (a), ... a,»)*, and a, = On 4m-1/%,-12 we have 





6S ial ( ie Mae 1 
= ae ————_-— eee ee no eo an ae 
” a, a ee akak., ... ak 








n n+l n“n+1 n+m-1 
1 1 i 
+ K + K P28 oe k is 
Oy Ay Oe Tan +] Oy Ty 4) n+ m-1 
+ ( ; + ; + + 
an afaar,, 
=. (- + falas i + opti v : — : + ) 
= n ee k k eee k k k = a eee 
G.7Mn Gage; ee Oy. a 


a rational number. 


Note that for the ordinary Engel series the condition a,,, 2 a, implies that for 
some ” sufficiently large a,,; = a, for alli 21. 


4. k*'® Power Series Related to the Liroth Series 


We could at this stage investigate expansions for real numbers whose terms 
take the form of the terms of the Litiroth series raised to a power. However, we 
consider instead a similar type of algorithm which leads to an expansion of 
simpler form, yet where the digits satisfy similar conditions. In particular, 
by setting kK = 2 in the results below, we obtain a series: expansion for real 
numbers with the appearance of a simplified Liiroth series. 


Theorem 4.1: Let k > 0. Every real number A has a representation 
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1 1 1 
A=a,+—* ; 


a aka (a,a,)* ” ( yk 
l 1%2 1727 43 214543) ay 





where: 
if ki> 1, then a, 2 2, a; 2 1 for +21, and a; 2 2 infinitely often, 


if 0 <k< 1, then Q34) 24; 2 2 fori 21. 


Proof: We derive this result from the following algorithm. Given any real num- 
ber A, let A) =A-a,0O<A, <1. Then we recursively define 


1 
Qa, = 1+ [=] for nv = 15. A, > 6, 
n 
where 
Ans) * aka, = qk-} for A, > 0. 


Applying this algorithm repeatedly, we obtain 


l Ay 
A= ay + A, = Fy FS eM ans 
y 
1 1 1 An+] 
ig 0 EEK 6 oes ee me Gn, 
ay aa, (a, sais a,-1) a, (a, eee Q,) 


Now a, = 1 + [1/A,,] implies that for A, > 0, 


A, > me 
n é 
and provided a, > 2 
A, < : 
me 


Now 0 < A; < 1 implies that a, 2 2 and Ap = ata, = ay > 0; thus aj 2 1. Con- 


tinuing this process inductively, we see that A, > 0; hence, a, 21 for all n. 
Consider the case k > 1. Suppose now that a, = 2; then 


ak _ a 
a =~ 1 a, - 1° 
Now if A,4) < 1, then a,4; 2 2. Otherwise, a,,; = 1 and A,42 = An+)} - 1. Con- 
tinuing this process, we see that after at most [ak-1/(a,, - 1)] steps with 
be A =A -1=A - 1, 


we must reach a stage at which A,,; < 1 and a,4; 2 2. Hence, there exists a se- 
quence of integers m, = 0 < nm, < n3 < ... such that 


anti = n+it+l n+7 n+l 


0 < Ay, +1 S 23 Ay +) 22, and a, = | 


for all other n> 1. Then 


An. 4) ' 1 


0 < —_— | 


.* Ba-n 
(aa, +++ ay.) 2k@-)) 


and so S,, + A as <t + », where 
1 1 Ans 
a = ao +— Fees ¢ ee = Sareea 
ay (a, +++ 4,1) ay (a, +++ Gy) 





1989] 57 





REPRESENTATIONS FOR REAL NUMBERS VIA k‘® POWERS OF INTEGERS 


Now let m;_, <n <n;. Then 
Sy, <S, < S,., andn +o iff 1 + », 
t-1 i 


So S,>A as n > », and the series converges. For the case 0 < Ks 1 te a, = 2 
then 


Ag eS 





+1 ‘iz - 
’ Ans) ak : 

Since k < l, qh k < dy and: as a, 2 2 we deduce that as a, 2 2 tor abi.in 2 Fs 
Thus, 


A 
me < as >Oasn>+o, 


(a, ... a,)* ake 


and again the series converges to A. 

We note that by setting k = 2 in Theorem 4.1 we obtain the expansion 

ae ee * 5 
A ATA (ag) “ag 

where a, = 2, a; 2 1 for 7 2 1, and a; 2 2 infinitely often. In many ways this 
could be regarded as a simplified version of the Liiroth series. In addition, 
we shall show shortly that, as in the Liiroth case, A is rational if and only if 
the digits in the expansion become periodic. 

A second algorithm for k 2 1 leads to a "binary" series of this type where 
the digits a; are equal to 1 or 2, for 7 2 1. 


Theorem 4.2: Let k > 0. Every real number A has a representation 


7 1 1 1 
a, +s 


k k 
ay a}a, (a,a5) a3 





where: 


it K-> 1,. then i =-2, 148, 42 for ~ 2 2, and a; = 2 infinitely often, 
1 


if 0<k< 1, then a, = 2, <6 gal-* for ¢ 2 As 

and a; 2 2 infinitely often. 
Proof: We use the algorithm of Theorem 4.1 except that now we let a, + A- aq» 
1 < A, < 2, and 


2 
a, = 1+ va for n>1, A, > 0. 
n 


In the same way as before, we can show A, > 0 anda, 2 1 for all n2l. 
Also in this case 


2 

Aes 

an 
which implies 

A 


eee eK x pk=l k-1 
n+l anA, ay, = an . 


It follows that for k > 1 


a =l1+ : 
n+1— 
A, 








< 1+ a 
k-1 
+1 ay 
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Thus, 1 < a, < 2 for m2 2 and, as 1 < A, < 2, a, = 2. Also, provided a, = 2 
(the case qa, > 2 cannot occur from the above), we get 
2 
A Se 2s 
a, =! 
and 
2an k-1 k-1 
An+] < <4 a = 3.2 > 
since we assumed a, = 2. Now in the same way as in the previous theorem, after 
at most [3.2*-!] steps of An+i = 1, we reach a stage at which A,,; < 2 and ay4; 


= 2. The convergence of the series for k 2 1 is now shown in exactly the same 
way as in the previous theorem. However here, unlike that case, the sequence 
{A,} is bounded above for all nm by a fixed constant as well. The proof for 0 < 
k < 1 is the same except that we now have, for n > l, 


44 2 3 x it FO. 


We consider now the expansion of rational numbers via these algorithms when 
k is a positive integer. We show that as in the previous section A is rational 
if and only if A has an expansion in which the digits become periodic. 


Proposition 4.3: Let k = 1, 2, 3, ... . The digits in the k*® power expansions 
of Theorem 4.1 (or Theorem 4.2) become periodic if and only if A is rational. 


Proof: First suppose that the expansion is periodic, that is, eventually 


ay a an+m- 


Then with the notation 





. 1 l " 1 
Be ae ee iB a GO ( " 
ay a\25 ay eee an-2 an-} 


i and a, = (Qa «-. ar, a, = Dean j/%_ ns we have 














1 1 1 1 
A=X, + — + ——— + oes + — 
’ : ( Ka 
On-1 (En AnIn+] Gy 2° 2n4m-2 n+m-1] 
1 1 1 
+ + ae ee + eee + esha R NET eae 
MA, %A,a,4) Wyld, 46% Ogee @ cei 
+ : ; + 
2 Zak a 
O5Qn A40,d,4) \ 
XxX, + : ( > ; — 
= n —" k eee 
On-1\%n Antn+] 
: (1 + : - Dos 
r ( a Ka i a 
hy +06 -Wgeeal Genns * 


a rational. 


Conversely, suppose A, = p/q is rational (with p, q © N). Then, since KEN, 
each A, is also rational, with 


aia eS a 
Ay * @ Aly ~ Guy" G eh. eF}) a 
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where DE Z, a, p, © N. Now in the case of the second algorithm (Theorem 4.2) 
“2 k=-1 
0<A,=p,/q < 3.2 : 


Thus, every 


.-s 3.2k-1qg 
es es 


and we deduce that the expansion becomes periodic in the same way as in 
Proposition 3.3. In the case of the algorithm of Theorem 4.1, we do not have 
a fixed bound for A,. However, when A is rational, 





1 q 
mie) eo) ee sg 4k 
ay La pe q 
as p, 2 1 for A, > 0. Using the fact that any A, for which a, = 1 is bounded 
above by an A,, for which a, 2 2, and that (for a, 2 2) 
qk-1 
An <—— 
at 
it follows: that, for all n 21, 
A, < (q+ 1)k7}. 
Thus, every 
{3 2 (q + a 
iE 45S) a eee 
q q q 


and again we can deduce that the expansion must eventually become periodic. 

In summary, we have found new classes of representations for real numbers 
that are related to the classical series of Sylvester, Engel, and Liiroth. In 
many cases, the expansions require very mild growth conditions on the digits 
and share with the Liiroth series the property of begin periodic when a number 
is rational. Unlike the classical series, however, the expansions for real 
numbers with kK # 1 are not unique, slightly different algorithms yielding 
series with the same properties, but with different digits for the same real 
number A. 
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1. Introduction 


In Chapter 3 of his second notebook [1l, p. 164], Ramanujan defined polyno- 
mials 
F=% 
A,(x) = Yia(r, k)x2P-k-l (p > 2) (1.1) 
k=0 


with A,(x) = x. The numbers a(r, k) are integers such that a(2, 0) = 1 and, 
for r= 2, 


a(r +1, k) = (rv - loha(v, k - 1) + (2r - k - l)a(r, k). (4.2) 


Also, a(r, k) = 0 when k < 0 or K > vr - 2. Properties of A,(x), and the moti- 
vation for defining them, are discussed in [1l, pp. 163-166]. Included in that 
reference is a list of the polynomials A,(x), 1 < r < 7, and the following theo- 
rem: 

r-2 

dar, k) = A,(1) = ( - 19771. (1.3) 

k=0 


In [3] it was shown how a(r, r - kK) can be expressed in terms of Stirling 
numbers of the first kind, and the following special cases were worked out: 


a(r, 0) = 1*3*5 eee+e (2r — 3), (1.4) 
a(r, 1) = [1*°3*5¢ «++ © (2r - 3)](r - 2)/3,; (1.5) 
a(r, r- 2) = (vr - 2)! (1.6) 


We note here that it is easy to prove by induction that 
a(r, 2) = (rv - 3)(r - 2)(r - 1)5° 7 e+ © (2r - 5)/3. 


The main purpose of the present paper is to prove congruences for a(r, k) 
(mod p), where p is a prime number. As an application of some of these congru- 
ences we prove Ap (x) /a?*} and Ap-,(x)/x?P are irreducible over the _ rational 
field. We also determine, for all r, the least residues of a(r, k) (mod 2), 
(mod 3), and (mod 4). For each r we find the largest k such that a(r, k) # 0 
(mod p), and we make a conjecture, based on computer evidence, about the 
smallest k such that a(r, k) #0 (mod p). We also conjecture the following 
periodicity property: 


a(r + (p - 1)p, k + (p - 2)p)=a(r, k) (mod p). 


This has been verified for all primes p < 251. A few other results and conjec- 
tures are given for moduli not necessarily prime. 


a Congruences (Mod P) 


Theorem 2.1: For any prime number p, 
a(p, k) = 0 (mod p) (k= 0, 1, .-+3 p - 3); 
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a(p, p - 2) = 1 (mod p). 
Proof: In [l, p. 164] we have 

A,(a) = x(r - 2)A,_1 (2) + 2 (” Nay ()Ay 4 (2) » (2.1) 
and hence | 

Any,(") = -@Ap(x) + @Ap(x)A, (x) = (x? - x)Ap(x) (mod p). 
Comparing coefficients of 2p -k+1, we have 


a(p + 1, k) = a(p, k) - a(p, kK - 1) (mod p). (2.2) 


From (1.2) we have 


a(p + 1, k) = -a(p, k - 1) - (kK + 1)a(p, k) (mod p). (2.3) 
Combining (2.2) and (2.3), we see that 
(k + 2)ates ®) = 0 Cand p) (k= 0, ssu5 p= 2). (2.4) 


The theorem now follows from (2.4) and (1.6). We note that results similar to 
Theorem 2.1 have been proved for the Stirling numbers [2, pp. 218-219]. 


Theorem 2.2: For any odd prime number p, 
Bip--.i, KX) = 0 Ged pp) (K=O. 2s .ses BD = 4), 
ap = 1, p = 3) = (Pp = 3)! “@iedp). 
Proof: From (1.2) we have 
a(p, k) = -2a(p - 1, k - 1) - (kK + 3)a(p - 1, k) (mod p). 
Thus, by Theorem 2.1, 
(k + 3)a(p - 1, k) = -2a(p - 1, k= 1) (mod p) (kK = 1, ..., p = 3). 
Since (225) 
a(p - 1, 0) = L« 3+ «++ «(2p < 5) = 0 Guet p) for-p > 3, 
the theorem follows from (2.5) and (1.6). 


Theorem 2.3: For any odd prime number p, the polynomials 
Ap(x)/2?*! and Ayu, (a) /x? 


are irreducible over the rational field. 


Proof: Assume p > 2. We know 
dpe kX) 28 (ued p) (eh = 6, |, scene — Bs 
a(p, 0) = 1+ 3+ +++ ©(2p - 3) #0 (mod p?), 
a(p, p- 2) = 120 (mod p). 


Thus, A p(x) /xP*! is irreducible by Eisenstein's Criteria. The proof is similar 
for A, (x) /x?. 


We note here that Theorem 2.1 could be generalized by using p/, j 2 1, in- 
stead of p. Replacing p by pJ in the proof, we have 


a(pi, k) = 0 (mod p) (k # -2 (mod p)). 
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Theorem 2.4: If p is an odd prime and if m 2 p, then 
a(m, k) = 0 (mod p) (kK =0, 1, ..., p - 3). 
If m > p, then 
a(m, p- 2) =a(p +t, p - 2) = le 3 «++ e(2t - 1) (mod p). 
Proof: We use induction on m. The first part of the theorem is true for m = 


Assume it is true form=p, p+1, ..., yr. Then, by (1.2) we have, for kh 
yaa aa 


a(r +1, k) = 0 (mod p); 


therefore, the first part of the theorem is true for all m > p- Now, by (1.2), 
if G > 0, then 


a(p + t; p = 2) = (2¢ = l)a(p + t - 1, p — 2) 
= le Fes 8 (26>= Vatpsip —-2) 
S le Je eee (24 - 1) (mod p). 
This completes the proof. 
We note that when t 2 (p + 1)/2, 
alp + tsp = 2) = O Ged p) (p> 2). 


We also note the following special cases of Theorem 2.4. For k 
Pe in Ne 


a(p + 1, Kk) = 0 (mod p); a(p + 2, k) = 0 (mod p): 

a(p + 1, p - 2) = 1 (mod p); a(p + 2, p - 2) = 3 (mod p); 

a(p +1, p - 1) = -1 (mod p); a(p + 2, p - 1) = -2 (mod p); 
a(p + 2, p) = 0 (mod p). 


Theorem 2.5: Let p be an odd prime. Then, for k = 0, l, ..., 2p - 5: 
a(2p, k) = 0 (mod p); (2p - 1, k) = 0 (mod p); 
G(2p, 2p - 4) = 1 (mod p); a(2p = ly 2p ~ 4) = 1 (mod) p); 
at2p;. 2p. = 3) = —2 rod py; a(2p - 1, 2p - 3) = 0 (mod p); 
a(2p; 2p = 2) = 0 (mod p). 


Proof: We know by (1.6) and Theorem 2.4 that 
a(2p, 2p - 2) = 0 = a(2p, p - 2) (mod p). 


From (2.1) we have 


Aon 4, (@) = (-w + @7)A,, (a2) + 2@Ap(x)A,,\(e) (mod p). 


Thus, by Theorem 2.1 and Theorem 2.4 (with m =p +1), 

Ay, 4;t@) = Ge + x*}A.., (2) + 2x¢2P*5 — 242P+4 (mod p). 
thaeiennce (2.6) gives, for k + 2p - 3, 2p - 4, 

a(2p + 1, k) = a(2p, k) - a(2p, k - 1) (mod p), 
and from (1.2) we have 

a(2p +1, k) = -(k + 1l)a(2p, k) - a(2p, k - 1) (mod p). 
Combining (2.7) and (2.8), we have, for k # 2p - 3, 2p - 4, 
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(k + 2)a(2p, k) = 0 (mod p). 

For k = 2p - 3 and k = 2p - 4, (2.6) and (2.8) give 
(2p - l)a(2p, 2p - 3) = 2 (mod p), 
(2p - 2)a(2p, 2p - 4) = -2 (mod p), 


and we see that the congruences for a(2p, k) in Theorem 2.5 are valid. Now, by 
(1.2) and (1.4), we have 


a(2p, k) = -2a(2p - 1, k - 1) - (Kk + 3)a(2p - 1, k) (mod p), 

a(2p - 1, 0) = O (mod p). (2.9) 
Thus, a(2p - 1, k) = 0 (mod p) (kK =0, 1, ..-5 p - 4), 
and by Theorem 2.4, 

a(2p - 1, p - 3) = 0 (mod p). 

It is now clear that the congruences for a(2p - 1, k) follow from the congruen- 
ces for a(2p, k) and (2.9). This completes the proof. 
Theorem 2.6: If p is prime and m 2 2p, then 

a(m, k) = 0 (mod p) (k=0, 1, ..-, 2p - 5); 

a(m, 2p - 4) = a(2p + t, 2p - 4) = 1*¢3¢ +++ (2 +1) (mod p). 
Proof: We use induction on m. The theorem is true for m= 2p. Assume it is 
true for m= 2p, 2p +1, ..-, yr. Then, by (1.2), we have 

a(r +1, k) = 0 (mod p) (k= 0, 1, ..+s 2p = 5)$3 


therefore, the first part of the theorem is true for all m2 2p. By (1.2), we 
have, for t > 0, 


a(2p i Gs 2p = 4) 


(2t + l)a(2p + t = 1, 2p = 4) 
3°5¢ +++ ©(2t + 1)a(2p, 2p - 4) 
Le 3*5¢ +++ ©(2t + 1) (mod p). 


tt 


This completes the proof. 


Using the same sort of proof as the proof of the first part of Theorem 2.5, 
we can show, for p > 2, 


a(3p, k) = 0 (mod p) (kK =0, 1, ..., 3p - 73 k # 2p — 2). 
The case a(3p, 2p - 2) has not been resolved. We indicate with Conjecture 1 in 
Section 4 what the general situation appears to be. The next theorem deals 
with a related problem, namely, the problem of finding the largest k such that 
a(r, k) # 0 (mod p). 
Define g(p, r) to be the largest k such that a(r, k) # 0 (mod p). 
Theorem 2.7: Let r be a positive integer, yr >2. Write 
r=2+(s(p-1) +t)ptu 
withs 20,0<t<p-2, andQO <susp-1. Then 
g(p» r) = m= sp(p - 2) + t(p- 1) +u. 
Furthermore, 


a(r, m) =ul (p- 2)!/(p - 2 - t)! (mod p). 
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Proof: We give a brief outline of the proof by induction on s, t, and u. Note 
that showing a(r, k) = 0 (mod p) for all k > mis simple, and we omit the de- 
tails. The recurrence relation (1.2) is the main tool in all of the following. 
The theorem is certainly true for s = ¢t=xu=0. For fixed s and t, induction 
on u is straightforward. Note that the induction applies to arbitrarily large 
u; the statement of the theorem restricts u to the nonzero values of a(r, m). 
If the theorem is true for some fixed value of s, u = p - 1, and some value of 
t, then it is not hard to show that the theorem must be true for the same s, u 
= 0, and the successor of t. By induction, if this theorem is true for some s 
and for t = u = 0, then it is also true for that s and allO < t < p - 2 and 
OS eS pm Bs 


Now suppose the theorem is true for some s and all t and u such that 0 < 
tsp-Zaned0susp-=- 1. Let 


ry = 2 + (s(p - 1) + (p - 2))p 
and let 


My = sp(p ~ 2) + (p - 2)(p - 1). 
Then, putting ¢ = p - 2 in the induction hypothesis, we have 
a(ry + Us, my + u) = u! (mod p) for O<usp-l. 
Since ry - 2 = 0 (mod p) and 2rp - (my + 1) - 3 = O (mod p), we must have 
a(ry> My - 1) = 020-0! (mod p). 
Now induct on u to show that 
a(rg + us, m t+u- 1) = ueru!t (mod p). 
Finally, we can conclude that: 
a(Yy) + Pp» my) + p) = p! = 0 (mod p); 
a(rog +p» m +p- 1) =p-p!=O0 (mod p); 
a(2 + (s + 1)(p - 1)p, (s + 1)p(p - 2)) = a(rg + Ds my + p - 2) 
= 0+alry +p - 1s ag +P - D+ 4-0 - 3)* ay +p - 1. my + P= 2) 
(p - 1) °(p - 1)! = 1 (mod p). 


It follows that the theorem is true for t = u = 0 ands+1. By induction, the 
theorem is true for alls 2>0,0<t<p-2, andOsusp-l. 


The proof of the following theorem follows the same lines as the proof of 
Theorem 2.2. 


Theorem 2.8: Let p be prime, p > 3. Then, for 1 < t < (p - 3)/2; 
a(p - t, k) = 0 (mod p) (k=0, 1, .--5 p - 2t - 2). 
Ror 2°S-t-'s p= "is 
a(2p - t, k) = O (mod p) (k =0, 1, ..-, 2p - 2¢ - 2). 
For example, using Theorem 2.8, Theorem 2.2, and (1.2), we have, for p > 5; 
a(p - 2, k) = 0 (mod p) (kK =0, 1, ..+5 Dp - 6)5 
a(p - 2; p - 5) = -(p - 4)!/3 (mod p), 
a(p - 2, p - 4) = (p - 4)! (mod p). 


3. Congruences (Mod 2), (Mod 3), and (Mod 4) 


In this section we first determine when a(r, k) is even and when it is odd. 
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Theorem 3.1: 
Gt? ,.70) = 1. (nod 2) ((@°S 2); 
aise, 1) =» Gaed 2). <e-> 3). 
Gtr, kK) = 0 \Gnod 2) (k > 1). 
Proof: The congruences for a(r, 0) and a(r, 1) are clear from (1.4) and (1.5). 
By (1.2) we have, for k > 1, 
a(2r, k) = (kK + l)a(2r - 1, k) (mod 2). 
If k is odd, we clearly have a(2r, k) is even. .If k is even, then 
a(2r, k) = a(2r - 1, k) (mod 2). 
And by (1.2), since k - 1 is odd, 
a(2r - 1, k) = a(2r - 2, k) (mod 2). 
Thus, 


a(2r, k) = a(2r - 2, k) = +++ = a(kK + 2, k) = k! = O (mod 2). 
Now since 
a(2r + 1, k) = a(2r, k - 1) + (k + L)a(2r, k), 
we have a(2r + 1, k) is even if k > 1. This completes the proof. 
The patterns (mod 4) and (mod 8) are suggested by the computer and can be 


proved by induction on r. For (mod 4) we have the following congruences. 


Theorem 3.2: a(r, k) = 0 (mod 4) for all k except: 


_ ipl Gnod 4) 2f-r 
eee ae iG (mod 4) if r 


1 ‘or .2  ‘Gnod’ 4), 
0 er 3 (mod 4), 


ler 3 Gned 4), 


a(r, 1) = 0 teed ‘&), 


_ (mod 4) if r 
2 (mod 4) if r 


Wt ou 


a(r, 2) = 2 (mod 4) if r = 0 (mod 4), 
a(r, 3) = 2 (mod 4) if r = 1 (mod 4). 


Theorems 3.1 and 3.2 suggest the following, which can be proved by means of 
(1.2) and induction on n. 


Theorem 3.3: If k = 2n, then a(r,; k) = 0 (mod 2”). 


To prove congruences (mod 3) we need the following lemma, which is a spe- 
cial case of Conjecture 4 of Section 4. 


Lemma: For r > 2, a(r, k) = a(r + 6, k+ 3) (mod 3). 


Proof: The lemma is true for r = 2, since 
ats, 3) = 1. God) 3), 
a(8, k) = 0 (mod 3) if k # 3. 
Assume it is true for r=m- 1. Then, by (1.2), 


a(m + 6, kK + 3) = (m - 2)a((m - 1) + 6, k + 2) 
+ (2(m - 1) - 1 - k)a(m - 1, k) 


(continued) 
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= (m - 2)a(m- 1, kK - 1) + (2(m - 1) - 1 - k)a(m - 1, k) 
= a(m, k) (mod 3). 
Theorem 3.4: a(r, k) = 0 (mod 3) for all k except: 
ale, Wotiod 
: 2 


a(r jes4}) = (ry +1) (mod 3) if r # 0 (mod 6), 
. 2 1 (mod 3) if r = 0 (mod 6). 


1 (mod as PZ 25 


Proof: Suppose r = 2 (mod 6), i.e., r = 63 + 2. Then, by the lemma, 


a(r, k) = a(6(j - 1) + 2, k - 3) = +++ = a(2, k - 37) (mod 3). 
Thus, 


a(r, k) 


Wi 


{9 (mod 3) if k + 33, 
1 (mod 3) if k = 37 = (r - 2)/2. 


The other cases of r (mod 6) are handled in exactly the same way. 


4. Conjectures 


Theorem 2.4, Theorem 2.5, and information given by the computer suggest the 
following conjectures. 
Conjecture 1: For all integers t and positive integers h such that h +t 21, 
a(hp + t, k) = 0 (mod p), k = 0, 1, ...5 A(p - 2) - I, 
a(hp + t, h(@p - 2) = 1*3¢ -++ «(2¢ + 2h - 3) (mod p). 


For t => 0, Conjecture 1 has already been proved in Section 2 of this paper 
for h=1, h = 2. If we try induction and assume true for h =m - 1, we can 
show, as in Theorem 2.1 and Theorem 2.5, 


(k + 2)a(mp, k) = 0 (mod p) (k = 0, ..-s m(p - 2) - 2). 
Thus, the proof depends on showing 

a(mp, k) = O (mod p) if k = -2 (mod p). 
The rest of the proof, for t > 0, would then follow. The cases t 2 0 have been 
verified by the computer for all primes less than or equal to 251. The case 


h +t = 0 leads to the next conjecture. 


Conjecture 2: Let p be any prime. 
(i) Let h be any nonnegative integer. Then 
a(2 + hp(p - 1), m) = O (mod p) if m = hp(p - 2). 
(ii) Let h be a nonnegative integer, h # 0 (mod p). Then 
a(l1 + h(p - 1), m) = O (mod p) if m # h(p - 2). 
(iii) Let h be a nonnegative integer, h # 0 or p - 1 (mod p). Then 
a(h(p - 1), m) = O (mod p) if m# h(p - 2) - 1. 
By Theorem 2.7, we know: 
(i) a(2 + hp(p - 1), hp(p - 2)) = 1 (mod p). 
(ii) Let h=sp+t, 1<t<p-1, s 2 0, then 
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a(l + h(p - 1), h(p - 2)) 


a(2 + (s(p - 1) + (¢ - 1))p + (p - 1) - £); 
sp(p - 2) + (é - l)(p - 1) + @ - 1) - &) 


(p- 1-t)!*° (p - 2)1/( - 2 - (€- 1))! = 1 = (mod p). 
(iii) Let h=sp+t, 1<t<p- 2, 8s <0, then 
a(h(p - 1), h(p - 2) - 1) 


a(2 + (s(p - 1)+ (t - 1))p + (—p - 2) - £); 
sp(p - 2) + (t - 1)(p - 1) + (p - 2) - £) 


(p-2-t)!*(-2D!/(p-2- (€-1))!=@-2D/M-1-t). 


The authors are grateful to the referee for suggesting the next conjecture. 
Part of this conjecture would follow from Conjectures 1 and 2. 


Define f(p, r) to be the smallest k such that a(r, k) # 0 (mod p). 
Conjecture 3: Clearly f(p, r) = 0 if r<(p+1)/2. Thus, for r 22: 


(i) f(p, r) = (p- 2) 5 S 4 - lif (p- 1)|r; 





(ii) f(p, r) 





(p - 2 5 - 4 if r=t (mod p-1),1<t< (pt+1)/2; 


IV 





(iii). f(p, r) 2 (p- 25 2 4 + 2¢ if r=t+ (p+ 1)/2 (mod p - 1), 
1<t< (p - 5)/2. 


In some cases, f(p, r) is larger than the formula given in (iii) above. 
For example, f(1l, 17) = 13, f(1l, 48) = 42, f(13, 22) = 19, and f(13, 68) = 59 
are larger by 2, and f(41, 350) = 334, f(43, 1743) = 1703, f(61, 2152) = 2111, 
and f(67, 2038) = 2002 are larger by 4. It appears to be difficult to predict 
when f(p, vr) will be larger than the formula or by how much it will be bigger. 


There are many cases where f(p, r) is larger by 2 or 4, and we suspect the 
formula could be off by even more for very large primes. 


Conjecture 4: If p is any prime, then 
a(r + p(p - 1), m+ p(p - 2)) = a(r, m) (mod p) for any r = 2, m2 0. 
Because of the recursion formula, (1.2), it suffices to show that 
a(2 + p(p - 1), m+ p(p - 2)) = a(2, m) (mod p) 


for all integers m. In this manner, Conjecture 4 has been proved on the 
computer for any prime p < 25l. 


Conjecture 5: If p is any prime and m 2 0, then 


a(r + p"™(p - 1), m+ p™(p - 2)) = a(r, m) (mod p”) 
for all sufficiently large r. 


Conjecture 5 has been proved on the computer for p” up to 2, 5 Fah 
115, 13%, 172, 19%, 23%, 292, 312, 372, 417. 


Conjecture 6: If p is an odd prime, then 


(i) a(r, 0) = p (mod 2p) for r = (p + 3)/2; 
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O (mod 2p) if r is even 
p (mod 2p) if r is odd 


(iii) if m= 2 and r 2 (p + 3)/2, then 


(ii). ate, ty = { and r = (p + 3)/2; 


a(r + p(p - 1), m+ p(p - 2)) = alr, m) (mod 2p). 


Conjecture 6 can be proved to be true if Conjecture 4 is assumed to be 
true. Similar conjectures for other composite moduli also seem to hold, but 
are more complicated to state. 


as Concluding Remarks 


Apparently not much is known about the numbers a(r, k). It would be useful 
if a generating function and a combinatorial interpretation were found. Also, 
it appears difficult to find values of A,(x) for x = 0, x = 1. We remark that 
it is easy to find derivative formulas for A,(x), however. It follows from 
(1.2) and the definition of A,(x) that 


x3At (x) = An, (0) - (r - 1)x,(x), 
(9) 


and thus it is easy to find a general formula for 4A, 
by (1.3) and the above comments, 


Af(1) = r? - (ry - 1)” 
A™(1) = (xv + 1)[(r + 1)” - 2r? + (vr - 1)*). 


(7). For example, we have 
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Let {x,},>,] be a sequence of positive real numbers. The q-sequence corre- 
sponding to {zx,,} is defined as having as its n™ term 


_ ent] 
Gn+1 iP Ly, 





for all integers m => 1. One of the purposes of this note is to compare the se- 
quence {z,,} with its corresponding g-sequence {q,} so that conditions imposed 
on one of them will yield results concerning the other. 


1. Example: Consider the Fibonacci sequence {x,,} defined recursively by 
Sly Sy le ey * SF PEM 2 


It is well known that the q-sequence corresponding to {x,} converges to the 
real number (1 + ¥5)/2. This example shows that divergent sequences {x,} can 
have corresponding q-sequences that converge. On the other hand, examples can 
be found of convergent sequences of quotients of convergent sequences. See 
Theorem 5 below. 


Whenever a sequence {x,} is defined recursively, say 
“=a, %, =b, x, = f(x,» x,_,) for n 2 2 


and positive numbers a and b, let S(a, b, f) denote the corresponding q- 
sequence, where f is a nonnegative function of two real variables which is 
defined and positive in the first quadrant and defined on the positive y-axis. 
If {q,,} converges, let z be its limit 
a= Lim dy. 

2. Theorem: Let {q,} = S(a, b, f) be the q-sequence corresponding to a sequence 
{x,}. If {q,} converges and f is continuous and positively homogeneous of de- 
gree 1 [f(Ax, dry) = Af(x, y) for »X > O], then zg satisfies the equation 


w2= f(w, 1). 


Proof: Since f is positively homogeneous of degree 1, it follows that 





xy 
Gn+19n oe , 7 = Itd.« 1). 


= 
Consequently, 2% = f(z, 1) must hold because of the continuity of f. 


3. Examples: (1) For the Fibonacci sequence, one has 
f(x,> x, 1) 7 F_ F Byis 
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and the limit 2 = (1 + ¥5)/2 satisfies f(z, 1) = ((1 + ¥5)/2)2, in agreement 
with the theorem. 


(2) Consider {q,} = S(1, 2, x, + 2x,_|). To find 2, one might 
want first to solve the quadratic equation 22 = z+ 2, whose positive solution 


is 2 = 2. Unfortunately, Theorem 2 as stated does not guarantee that lim 
2. If the limit exists, then _ 





= e. - 2e x 
_ See ae w=) _ E a ; 
Gey +g oe =1+ 2) = 1+ 2/4, 
implies 2 = 2. A procedure for finding the limit is presented in the next 
result. 


4. Theorem: Let b > 0 and c 20 be real numbers. If f(x, y) = be + cy, let 
{x,} be the sequence defined recursively by 


x, = p> Oi, Lo =q> QO; and a5 * Fy x, -1) for m = 2. 
Then the g-sequence S(p, q, f) converges to 
2 = (b + vp2 + 4e)/2 


independent of the initial values p and g. Moreover, the sequence {q,4) - q,} 


is either the constant sequence {0} or oscillates between positive and negative 
values. 


and hence, for nv 2 3, 


Bg Gn-1 si q 


Consequently, {q,,,-4,} is either the sequence {0} or oscillates between posi- 
tive and negative values. Also, 





Qn = b+ for n 2 3 


al 
implies that q,q,_, = bq,-, + ¢ > b2 + ce for n = 4, and hence 


i Wn-1 ~ In 
WnTn -) 


l€na1 ~ Inl S———(4,.) - Gal - 








If d = c/(b2 +c), then 0 < d < 1 and 


ldnsa - Mel * 2 * la, - G3] for » 2 3. 
Since 


eo 


va lay - q3\| 


n= 


converges, it follows that {q,} is a Cauchy sequence, thus it converges to some 
number z > 0. Theorem 2 shows that 2% = bz + c must hold; therefore, 


a Vp2 + 4e 
ee ee 


z 


O 
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One can, in some cases, compare the behavior of a given sequence {x,} with 
that of its corresponding g-sequence. It is to be pointed out here that the 


sequences referred to in the following result are not necessarily generated by 
recursion 


5. Theorem: Let {x,} be a sequence of positive numbers and let {q,,} denote its 
corresponding q-sequence. Then 

(1) If {q,} € e 9, then {z,} € a); hence {x,} € Co: 
(20), Co» and ec denote the spaces of summable, convergent to zero, and conver- 
gent sequences, respectively.) 


6. Example: Consider the sequence defined by 


a 
nan? - 1 
Since 4n* - 1 > n*, then {y,} € 2“). Define a sequence {x,} by 
t 


Tq + Yy) 


k2n 


n 


where t > 0 is a real parameter. This sequence is well defined, because the 
infinite product 


fl (1 + y,) 


converges. It follows that {1 + y,} is a q-sequence. Indeed, 


vn+1 - 


n 1+ 
roth Yx) 


hi * Y,? 


x 


for n 2 1. A simple computation shows that 


Marty) =-2ssi23- *, 
k21 


where the product on the right was shown by John Wallis (1616-1703) to have the 
value 1/2. 


The next result is an attempt to answer a question suggested by the pre- 
vious example: What sequences are q-sequences? 


7. Theorem: Let {y,} be a sequence of positive terms in 2“). Then, there ex- 
ists a l-parameter family of sequences {z,,(t)} for t > 0 such that 


ror ari. 2 1. 


Proof: In order to have that x ,,/x, = 1+ y,, one must solve the infinite lin- 
ear system 
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x, (1 + ¥,) -%* 


xv, (1 + Yo) - £3 = 


Set x; = t, an arbitrary positive real number. Then x» = t(1 + y,)» 23 = t(1 + 
y,)(1 + yo), and by induction, 


Luz) 2 F MW d+ Y,)- 


ken 


Therefore, the q-sequence corresponding to {x,(t)} is given by 


4,0) + Marty) a3 
Se Ys 
T(t) + TW (ty) | 


k<n-] 


This establishes the result. 


In Theorem 2 above, the limit 2 of a convergent g-sequence corresponding to 
a recursively generated sequence was shown to satisfy a functional equation 
involving the generating function for the original sequence. This generating 
function was required to be positively homogeneous of degree 1, continuous and 
nonnegative in the first quadrant. According to Theorem 4, if the generating 
function is the restriction of a linear form, then the limit of the q-sequence 
can be explicitly calculated and does not depend on the initial two terms of 
the original sequence. 

The result that follows explores the nature of the functional equation by 
characterizing the class of functions to which Theorem 2 applies and provides 
examples to show that the independence of the limit 2 of a g-sequence with 
respect to the initial terms p and q of the original sequence, which was one of 
the conclusions obtained in Theorem 4, no longer holds in the general case. 


8. Theorem: A function f: [0, ~) x (0, ~) + [0, ~) is continuous, positive on 
(0, ~)*, and positively homogeneous of degree 1 if and only if there is a con- 
tinuous function y: (0, ~) > (0, ~) which is such that 


(i) ye) = FU. t) for all ¢ € (0,.~) 


and (t) 
(ii) lim ~ exists and is finite. 
t>oa 


Proof: Suppose f:{[0, ~) x (0, ~) + [0, ~) is continuous, positive on (0, co) 2 
and positively homogeneous of degree 1. Set 
yCéy) = 7 Gls. t) fer & € (0; >}. 


Then y is continuous and positive, and 


tim 1) ~ yim f(4, 1) = FO, 1) 


t>+oa t>o 


exists and is finite, due to the continuity of f and the fact that f is posi- 
tively homogeneous of degree l. 
Conversely, suppose that y: (0, ~) > (0, ~) is continuous and that 


lim y(t) 


t>o 


exists and is finite. 


lim y(t) 


t > 


a= 
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and define f on [0, ~) x (0, ~) by setting 


xy(y/x) if x = 0 


f(z, y) = 
ya if 2: = OQ. 


f is continuous for x > 0 and y >0O since is continuous on (0, ~), (x, y) > 
y/x is continuous for x > 0 and y > O and the projection (x, y) + x is continu- 
ous everywhere. 


If S, = fix, y) 24 > 0; y > OF, 5, * {(O, y): y > 0} and if y, is a fixed 
positive number, then 


lim f(x, y) = lim wy = lim yi} * ot 
(x,y) >(0, yo) (x, y) > (0, ¥9) - (xz, y)*(sy¥.) | ¥ 
(x, y)€ Sy (x, y)€ Sy (x, yJ)ES, ” =f, Yo)> 
and lim f(x, y) = lim oy = ypu = f(0, yo). 
(x,y) >(0, yo) ¥*Yo 
(z,y) € So 
Thus, 


lim Tas y) = f(O, Yo)> 
(x, y)> (0s ¥9) 


and f is continuous at (0, y)). It has now been demonstrated that f is contin- 
uous on [0, ~) x (0, ~). 
The function f is also positively homogeneous of degree 1. For, if A > QO, 


: : 
ihe, tas hig if x +0 . : xy(4) if x z= 0 


(ay)a if x =0 A (ya) if x =0. 


Therefore, f(Ax, dAy) = Af(x, y). U 


Since Res '* f(x,, vp) implies 


6.24 
Ine. ~ JX? =) Nan)’ 


the question of the convergence of q-sequences is equivalent to examining the 


convergence of sequences {q,,} generated by choosing q, > 0 and defining q, for 
na 3 by 


ie) 
qn ios 


for some positive continuous function y on (0, ~) having the property that 


lim y(t) 


tro 





exists and is finite. Using this fact, examples can be constructed quite 
easily. The following examples, which were constructed in this way, show that 
limits of q-sequences depend in general on the starting values x, and xy. 


9. Example: Let y(t) = 1/t%. Starting with qz = 1, it follows that q, = 1 for 


all-n and lim q, = 1. However, if q, = 2, then it is easy to show that 
t>+o 


ae 


for n 2 2, so that {q,} diverges. This shows convergence is dependent on the 
starting values. In this example, 


f(x, y) = xy(y/x) = x3/y?. 
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ON A CERTAIN SEQUENCE OF QUOTIENTS OF A SEQUENCE 


x, and x, could be taken as x, = x, = 1 and x, = 1, x) = 2, respectively. 


10. Example: Let y(t) = 1/t and let q, = r > 0 be arbitrarily chosen. Then 


1 
“(e)-» 
q3 qo 
Similarly, q, = r for all n= 4. Therefore, {q,} is the constant sequence 
{Fo Pe occa Be seeds 


which converges to r. In this example, it is seen that each positive real num- 
ber is the limit of some q-sequence for the same generating function. Here, 


f(x, y) = y(“) = =. 


Starting values may be taken as x, = 1 and x, =r > 0. 
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THE ALPHA AND THE OMEGA OF THE WYTHOFF PAIRS 


J. C. Turner 


University of Waikato, Hamilton, New Zealand 
(Submitted February 1987) 


1. Introduction 
The Wythoff number pairs have been much discussed in the literature on 


Fibonacci integers (see [1] for references up to 1978). And in [2] and [3] M. 
Bicknell-Johnson treats generalizations of Wythoff numbers which provide number 


triples with many interesting properties. In this paper we present three 
different ways to generate the Wythoff pairs, and, with some trepidation in 
view of the extent of the literature on them, claim that these are "new." We 


emphasize the notion of "generation" (in contrast to "giving a formula"), and 
introduce Fibonacci word patterns [3] as a tool to define n-tuple generating 
processes. 

A determinantal relation for the Wythoff pairs is described, which makes 
further use of the word-pattern tools. 

In the final section we show how similar methods can be used to generate 
and study sequences of integer triples. Three examples are given, and each is 
an attempt to generalize aspects of the Wythoff pairs-sequence. 

It is clear to us that these tools and methods hold much promise for devel- 
oping a general theory of sequences of integer n-tuples which have structures 
related to Fibonacci word patterns. 


2. Notation and Definitions 


The main tool to be used below is the Fibonacci word pattern, which we de- 
veloped in [3]. We shall also use an operation of merging two integer 
sequences, and its inverse; we shall use the terms addmerge and submerge for 
these two operations. 


Definitions and examples 


To keep the exposition brief and readable, we now give somewhat informal 
definitions of the operations and concepts we wish to use. The examples will 
make the intended operations perfectly clear. 


(i) Ftbonaect word patterns (W,WoW, ... W, ..-) 


A word pattern is a concatenation of a sequence of words Wy> Wo» W535 


-, Wy» «+. + The words are formed using characters from a given 
letter set such as {0, 1} or {a, b, e@}. The basic word pattern is 
obtained by repeatedly using the concatenation recurrence 


Wie2 = W,W, 4, with W, = A, Wy = B. 
We shall denote the resulting pattern by F(A, B). Then: 
PAs B) = As Bs ABs BABS ABBABs sae % 


(N.B. The commas on the right should be removed; they are inserted to 
show the boundaries of successive words in the pattern.) 
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Examples to be used below: 


a= (1,0) =1010010100100... (using A=1, B=0); 
w= (1, 01) =1011010110110... (using A=1, B=01). 


These two binary word patterns are, respectively, the alpha and the 
omega referred to in the title of this article. The w pattern is named 


after Wythoff, for reasons which will become abundantly clear as the 
paper develops. 


We shall also use the tribonacci word pattern (with W, =A, W,=B, W3=C) 


F(A, B, C) = A; B, C, ABC, BCABC, ... (W, = Ee. 


n+3 +42) 


Set-sequences 


In [4] we introduced the following construction (though with a slightly 
different notation). Let {S5,} be a sequence of sets, and let {a,} be a 


sequence of integers. The set-sequence is formed using the following 
recurrence 


Ss =S US 
n 


n+2 + an» 


n+l 


with 5,, So being any given sets. The + operation is to be carried out 
as indicated by 


{8,5 Sos coos Sze eee} +a = (8) +s So + Ay woes 8; HAs «e}. 


Addmerging and submerging 


A merging operation, and its inverse, should be clear from the following 
definition and example. 


Let S and 7 be any monotone increasing sequences. Then the addmerge of 
S and T (written S ~ 7) is obtained by taking the multi-union of the 
two sequences and sorting them into monotonic increasing order. By 
"multi-union" we mean that integer repetitions are to be allowed. 


Example 
lee 66 8s he Ge Bi Bean OT ow 18; 35 Be ike ccks Toe 
Ban? @ thy 2 Se Ss B72 Bee Se My cvs 


The inverse of addmerge is submerge, which we shall write S~T. This 
operation simply removes the sequence S from the sequence S (all elements 
of 7, that is, which happen to be in S). 


Sequence notations 


We shall use either Greek letters or underlined, lowercase Roman let- 
ters to denote sequences; and will use or omit subscripts on individual 
sequence elements as is appropriate. The following examples illustrate 
our notation, and will be needed below. 


n = {l, 2, 3, .--} the natural numbers 


n* = {0, 1, 2, ...} the natural numbers with zero 

f = {1, 1, 2, 3, 5, ...} the Fibonacci integers {F,} 
£° © C85. 2s 35 Ss vee} © (hygy ls Bom Us. 25 cee 

Pi -4Ey By Sy Bid gg lh Wi Os cx 
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w, = {l, 3, 4, 6, 8, ...} first members of Wythoff pairs, 
equals {[na]} where a = $(1 + Y5). 
wy = {2, 5, 7, 10, 13, ...} second members of Wythoff pairs, 
equals {[na2]}. 
-{(*)- ene a 
we a = CMe relics) ... Wythoff pair-sequence. 


(v) Binary word pattern representations 
Let B be a sequence of 0's and 1's, say, 


DW bey Bes Bye ssc Bs, 2 OO OO: 
And let 


S)» So» S3o +++» Sy» +++ be an integer sequence. Then, 


= the subsequence from s whose elements are in the 
positions where the 1's occur in B. 


ty |H 


S x 


For example, if B= 0, 1, 0, 1, 0, 1, 
n*B=2, 4, 6, 8, 


estoy ERCRs 


3. Three Ways of Generating the Wythoff Pairs 


We now use our word patterns and sequence notations to give three differ- 
ent ways to generate the Wythoff pairs. 


(i) Use of the omega sequence 
Recall from Section 2 that the binary word pattern omega is 
i =O OL) =) 1-0) 4 Oa Om 01 


It may be observed that the 1's occur in positions w; = 1, 3, 4, 6, 8; 

-; and the 0's occur in positions wo = 2, 5, 7, 10, 13, ... . Thus, w 
contains all the information needed for producing the Wythoff pairs. 
Using the notations of 2(iv) and 2(v), we can write: 


si ea 2 (2 * FO, OL) 
= W n * F(O, 10) 
Note that we can do certain algebraic operations with sequences and our 
new notation. 
Thus, for example: 
n=) ™ W =n * EAGL; OY) + Os. 10) 
where + is mod 2 addition of elements. 


In [6] two methods of proof are given to demonstrate that omega [i.e., 
F(1, 01)] does in fact generate the Wythoff sequences 


w, = {[na]} and Wy = {[na2]} 
as claimed. 
(ii) Use of Fibonacct set-sequences 


In 2(ii) above, we explained the recurrence for generating Fibonacci set- 
sequences, viz., 


S =z=S US +a 


n+2 n n+1 n° 
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Let S_, = {0} and Sy = {0}, and {a} 2:£° = lp 2s 3a Se acs < They 
S, = {1}, S, = {3}, S; = {4, 6}, etce.; 


and it soon becomes clear that 
Us, =,. 
n=1 we =! 


Similarly, if {a,} = £" = 2, 3, 5, 8, ... and the same S_,, Sy) are used, 
the infinite union 


n=1 
Proofs of these assertions are given in [6]. 


Use of Ftbonacet magic matrices 


In [5] we decided that square matrices all of whose elements were Fibo- 
nacci integers, whose diagonal, row, and column sums were Fibonacci 
integers and, moreover, whose powers also possessed these properties, 
deserved to be called magic. The spectral radius of these matrices is 
a, the golden mean. The simplest such matrix is 


0 1 1 1 
A= . Note that A2 = 
1 1 1 Z 
The characteristic polynomial of A is a2 - a - 1; and of A? is \% - 3h 


+ 1, which has maximum root \ = a2 with a = $(1 + 75). 


Many properties of A have been noted in the literature, but the follow- 
ing relationships with the Wythoff pairs may possibly be new. We give 
them without full proof. 


Proposition (generation of w, the Wythoff pairs sequence): 


(A) w= {tna + mt1(1)}s 


where J is the 2 x 2 identity matrix, n is the natural number se- 
quence, and 


— {m,} = 0, ne ky 2% aa ds 4, 4, 55 6; 6; qs 


(N.B. The generation of m by a Fibonacci word pattern is given 
below. ) 


(B) w= {1,4 + mA?1(t)}, where 


i in} = be Peds) Saclay os 35 4, 4; 4, a3 
Note that (B) follows from (A) since m+ r = n and A2=A+TI. 
Generation of mand r 


The sequence m is generated as follows. Take the Fibonacci word 
pattern 


FPQly 2) ba dy Bas 212» P2225 


We can use the elements of this sequence as frequencies, drawing 
elements from the sequence n = 0, 1, 2, 3, 4, ... with these fre- 
quencies. This gives a natural extension to the star operation 
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which we defined in 2(v), in connection with binary words. Thus, 
we get 


m=n** F(1, 2) = 0, 11, 2, 33, 44, 5, 66, etc. 


as required. With this very useful extended operation (which in- 
cludes the earlier one, if 0 frequency is interpreted as "leave 
out"), we see that the sequence r in proposition (B) for w is: 


$* by Be Bi 2 Be BG Ae HOM «es 9S ME, DD. 
Corollary: n* * F(1l, 2) +n * F(2, 3) = n, since m+ r =n. 


The attractiveness of the method of generation of the Wythoff pairs 
just given lies in comparisons that can be made with the classical 
generation of the individual sequences. To spell these out, we 
note that w; = 1, 3, 4, 6, 8, ... is generated by [ma], and wo = 
2, 5, 7, 10, 13, ... is generated by [na2], where a is the golden 
mean. By comparison, the generation formula given in (A) above for 


» = (2)(5)z)(o)(r3) --: 


uses only mA, where A is a matrix having a as spectral radius, and 
a sequence n* * F(1, 2): and the sequence F(l, 2) has the same 
pattern as that other "a," the basic Fibonacci word pattern 
referred to in the title of our paper. 


4. A Determinantal Relation for the Wythoff Pairs 


The following interesting relationship is reminiscent of the well-known 
Fibonacci relation 


he Suut 
oe Fn+2 


Consider the Wythoff pair-sequence 
» = (2), 
an 22 
To simplify the notation, we write w, = a = {a;} and w, = b = {b;}. Then, 


oy = (2)(32)(33) = (RVG) : 


Let us pass a 2 x 2 window along this sequence and compute determinants as we 


go. Thus, 
t. -l, 2. -2, -2Z, ae -3, &, -4, -4, 


os 


a; i+ 


b; bes 


t 








6, -5;5 -5, os -6,; 9, -7; -7; EL. -8, 


There is clearly an interesting pattern to the sequence, but how can we capture 
it in a formula? It is here that our word pattern notation becomes really 
useful. Let us submerge the negative and the positive elements, to find: 


-l, -2, -2; -3, -4, -4, -5; -5; -6,; -7; -7; -8, 
and 


l, 3» 4, 6S. 8, a. Pas eee 
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Now we see that the negative sequence can be written as (-n) * F(1,2). And the 
positive sequence is just w). 


Thus, we state finally 


Proposition: 


(i) The determinants of successive Wythoff pairs are given by 


Aw = 








=w, ~ [(-n) * FCI, 2)], 





(with the addmerge ignoring minus signs). [J 


It might be said that to complete the above proposition we must give a 
precise formula for the m'® term of Aw, whereas we have given only a 
sequence generator. We shall do this later. As we said in the Intro- 
duction, we wish first to emphasize ways in which our notation can 
describe the generation of interesting sequences. Picking out par- 
ticular values of a sequence is always harder to do. In [3] and [6] are 
given many formulas for making that task easier, being results 
concerning counts of runs and runlengths of particular letters or 
integers in given Fibonacci word patterns. 


(ii) Formulas for the nt? term in the sequence of determinants are: 


w,);» when 7 € w, with n = [ia7], 
where a is the golden ratio; 
(Aw), = g& > 


-i, when [(7- 1)a2] <n < [ta2]. 


Proof: It may be seen that the positive terms in the sequence of 
determinants, namely, 


ty Be Su ha BA W. Be anes 


occur at positions 


: W, = 2, 5, 7, 10, 13, 15, BS}. ea, « 
This is fascinating in itself, and immediately explains the given 
formulas, because the positive terms occur when 7 = [ta2]. O 


[N.B. Because of the fact just noted in the proof, we could give the 
determinant sequence as 


dw = n * F(O, 10) ~ (-n) * F(1, 2). 
An immediate corollary of the fact that Aw includes the sequence 
(-n) * F(1, 2) 


is the following proposition on the representation of the natural 
numbers in terms of the Wythoff numbers. 


(iii) In terms of the Wythoff numbers, every integer V can be represented as 
follows: 


either uniquely as W = a,,,b; - b;,,4; using Wythoff pairs; 


t+1 
or in two ways using a run of three consecutive Wythoff pairs 
thus: NW = a,;,,b, - bp43%; = Oagodi4 ~ Pea2%ear- 


t+1 t+2 
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There any many other interesting things that could be said about the se- 
quence Aw. One more will have to suffice. Suppose we mark the sequence into 
words, each of which ends at a positive integer thus: 


(-1, 1) (=2; 2, 3) (-3, 4) (-4, -4, 6) (-5, -5, 8) (-6, 9) 
(-7, -7, 11) (-8, 12) etc. 


The lengths of these words have the pattern F(2, 3). And their totals follow 
the pattern 


0, =1., l; 2 -2, -2; 35 -3,; as -4, -4, 65 -5; -5, 3, a 
The first 0 indicates that the sum of the first two determinants is: 


If we consider the sequence of word totals, it appears that it will oscillate 
with increasing amplitude; and that the sum 


2 (Aw), 


t=1 


will equal zero infinitely often; but we have not established proofs of these 
observations. 


5. Generation of Other Pair-Sequences 


‘Any Fibonacci word-pattern which uses a binary letter-set can be used to 
generate a pair-sequence. For example, the alpha sequence 
a= F(1, 0) = 1, 0, 10, O10, 10010, OL0100I0, ... 


generates the following: 
eG -GG) | enone 


A question of interest now is whether a can be expressed in terms of the 


Wythoff pairs, and vice-versa. Using our word-pattern tools, we find as 
follows: 


g 
~ 
i 


bo 3-6°S EY 14 16 19 
= (1346869 11 12 14 16 17 19) ~ (69 12 17 ««s) 


= 0, ~ [w, * FCO, 01)) = w, * FU, 10); 


or 24S 79 42 19:19 27 .. 
—1757 0 OW. ~“ROR TT wd 
= w. ~ [w, * F(O, 01)]. 


Thus, we have: 
(“) ep wy [w, * F(O, ee) 
a = = 
= 7) Do oe [w, * F(O, 01)] 
By similar methods we can invert these equations thus: 


Oy Says [a, x F(O, an 


and 


Wo = ~ [ao * F(O, 10)] 
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And so the alpha pair-sequence can be expressed in terms of the omega 
(Wythoff) pair-sequence; and vice-versa. 


It is evident that by such means an infinite number of pair-sequences can 
be generated, and their properties studied by establishing relationships be- 
tween them and the fundamental Wythoff pairs. A new kind of number theory 
could be developed, based upon the sequences w, and w 2, and related to the 
“ordinary” number theory based on n through the functions [na] and [na2]. 


Finally, we give an indication of how these methods can be extended to 
study sequences of triples. 


6. Generation of Triple-Sequences 


We shall show, proceeding by examples and comments upon them, how to gen- 
erate triple-sequences in three different ways. The first uses Fibonacci word- 
pattern with letter-set {a, b, ec}; the second uses a tribonacci word pattern 
with letter-set {a, b, ec}; and the third uses a "magic" tribonacci matrix. 


(i) Use of a Fibonacci word pattern 
Consider the following word pattern: 
F(a, be) = a, be, abe, beabe, abebeabe, 
Listing the postttons of a, b, and c, respectively, produces the follow- 


ing triple-sequence: 


1 4 9 12 17 22 25 30 
2 a 7 10 13 15 18 20 
a 6 8 11 14 16 19 21 


lalo|p 
ul 


It will be noted that, as might be expected since the word pattern is 
Fibonacci, the three component sequences can each be expressed in terms 
of the Wythoff numbers. Thus: 


a W, * oO 
b = Wo » where a = F(l1, OQ). 


Then a, = (w)); >» b, = (w), » and c, = (Wo), + 1, where 
t= Cy.) ~ (@5),.; * I- 


Note also that a, b, c are each stiictly increasing sequences, they are 
non-intersecting, and their union equals Zt: all properties of the 
Wythoff pairs-sequence. Their proof is immediate from the way in which 
F(a, be) is expanded. 


Another interesting point is that the parity of the terms in a is alter- 
nately odd and even. And then, since the sum (b, + ¢@,) is always odd, 
we have the sum (a, + by + cy) also of alternating parity. 


The parity patterns, and more generally mod 3, mod 4, etc., patterns of 
elements of multi-sequences generated from word patterns would seem to 
be worthy of study. 
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Use of a tribonacci pattern 


Consider next the tribonacci expansion of F(a, b, ¢@) and the triple- 


sequence it generates through the positions of a, b, c in the resulting 
pattern. 


F(a, b, ec) =a; b, e, abe, beabe, cabebeabe, 


a 1 4 9 13 18 21 26 
Df pee 5 7 10 14 16 19 
ic 3 6 8 At 12 r5 17 


This triple-sequence again clearly has the property that each element 
sequence is monotone increasing, and the three sequences partition Z*. 


When we first studied this sequence, we hoped that we would find simple 
relationships between a, b, and c, respectively, 


{{nt]}, {[nt2]}, and {[nt3]}, 
where t is the positive root of the tribonacci equation 
x3 - 22-2 -1=0. 


This would have been an excellent generalization of the Wythoff pairs 
property whereby a = {[na]} and b = {[na2]}. Unfortunately, we have not 
been able to find such relationships, although there seems to be hope 
for relating Fibonacci word patterns to the sequences of first 
differences {A[nt] }, etc. To encourage the reader to search for these, 
we show the first few tribonacci triples (t = 1.839): 


[nt] 1 3 5 7 9 11 12 
[mt2] | = | 3 6 10 13 16 20 Oe PF bes 
[nt3] 6 12 18 24 31 37 43 


Use of a trtbonacet magic matrix 


Our third attempt to generalize the Wythoff pairs is to take a 3 x 3 
matrix which generalizes the "magic" properties of the 2 x 2 matrix used 
in Section 3(iii), and attempt to generate with it a unique sequence of 
triples whose members partition Z*. Once again we must confess that we 
have not found a fully satisfactory way of defining such a unique 
sequence; but in the spirit of the aims of this paper we believe it is 
worth presenting our attempt. 


The tribonacci magic matrix we shall use is: 


0 1 0 
Pee 30) ks 
1 1 1 
Note that the characteristic polynomial of T is -(A3 - \* - A - 1) so 


its spectral radius is t = 1.839. 


Note further that the row sums of powers of 7 are 


1 1 1 1 3 5 9 
T| 1 |; 72] 1 |, ..., which give | 1 3 5 9 IF he wand 
1 1 3 5 9 17 |{ 31 


with each element of the triples being in tribonacci sequence. These 
are the magic properties of TI. 
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The question we ask now is the following. In 3(iii)(B) we generated the 
Wythoff pairs using only the 2 x 2 matrices A and A? and coefficients 
from the sequences n * F(2, 3) and n* * F(1, 2). Can we generate a 
unique sequence of triples, a T-sequence, using only the 3 x 3 matrices 
T, 72, and T3, together with coefficients from sequences which can be 
defined in terms of Fibonacci word patterns? Furthermore, can we 
require the three member sequences of the triple sequence to be strictly 
increasing, and to partition Z*? If we can find such a T-sequence 
uniquely, it will constitute an excellent generalization of the Wythoff 
pairs sequence. 


Out attempt, down through the first twenty triples, is tabulated below, 
showing the triples horizontally for convenience. 







Triples (a, b, e) 
Elements 


Row sums of (pl + qT? + rT?) 
Coefficients 


P q 









OONNDUUNLSWNHeE CO 
WoOnnNnN DU LE DN & HE = 


W WN DDD FR RR RR RR RRR RK OO 


It will be noted that we have succeeded in advancing (a, b, c) thus far 
without increasing p, q, and r by more than 1 at each step. But, as we 
confessed above, we have not determined a formula for advancing the 
triple sequence indefinitely while satisfying all the requirements for 
generalizing the Wythoff pairs to triples. 


7. Summary 


In this paper we have defined word patterns, and various tools derived from 
them, in order to generate and study increasing sequences of irteger pairs and 
integer triples. 

A particular outcome of our study of pair-sequences as derived from Fibo- 
nacci binary word patterns was to show how all such sequences (and there is an 
infinite class of them) might be related to the Wythoff pairs. 
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It is hoped that we have convinced the reader that there is much scope for 
developing a number theory of integer pairs defined by binary sequence rep- 
resentations and using tools such as we have described. The title of our 
paper, namely, “The Alpha and the Omega of Wythoff Pairs," might suggest that 
all has now been said upon the pairs. In fact we claim the opposite—that this 
paper can mark a beginning of a broad development in their study and 
application. 

The path to a general study of triple sequences would seem to be a much 
harder (but nevertheless a most interesting) one to seek. 
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ELEMENTARY PROBLEMS AND SOLUTIONS 


Edited by 
A. P. Hillman 


Please send all communications regarding ELEMENTARY PROBLEMS AND 
SOLUTIONS to Dr. A. P. HILLMAN; 709 SOLANO DR., S.E.; ALBUQUERQUE, 
NM _ 87108. Each solution or problem should be on a separate sheet (or 
Sheets). Preference will be given to those typed with double spacing in the 
format used below. Solutions should be received within four months of the 
publication date. 


DEFINITIONS 
The Fibonacci numbers /F’, and the Lucas numbers L, satisfy 


PF = Ff 


n+l 


PF Pas OG fF, =} 
n 0 
and 


L = 1 hs La * 25 2,4. 


n+1 0 


PROBLEMS PROPOSED IN THIS ISSUE 
B-634 Proposed by P. L. Mana, Albuquerque, NM 
For how many integers n with 1 < n < 10® is 2" =n (mod 5)? 


B-635 Proposed by Mohammad K. Azarian, U. of Evansville, Evansville, IN 





For all positive integers n, prove that 


mel + > cat | < (mn + 2)"*!, 
k=1 


B-636 Proposed by Mohammad kK. Azarian, U. of Evansville, Evansville, IN 





Solve the difference equation 


Fes * (n + 1)x,, + r\(n + 1)3[n!(m! - 1)] 


n+ 


for x, in terms of i, x), and n. 
B-637 Proposed by John Turner, U. of Waikato, Hamilton, New Zealand 


Show that 
2 1 


Ze, Pes + aFy +1 


= 1, 


where qa is the golden mean (1 + V5)/2. 
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B-638 Proposed by Herta T. Freitag, Roanoke, VA 


Find s and t as function of k and n such that 
k 
me Sanmmeniae * 2y% 
i=1 


B-639 Proposed by Herta T. Freitag, Roanoke, VA 


Find s and t as function of k and m1 such that 
k 
he Ly -upsui-2 = 45 Z¢- < 
21 
SOLUTIONS 


No Fibonacci Pythagorean Triples 
B-610 Proposed by L. Kuipers, Serre, Switzerland 


Prove that there are no positive integers r, s, and ¢ such that (F,,F,, F,) 
is a Pythagorean triple (that is, such that Fe + F2 = Pee. 


Solution by Marjorie Bicknell-Johnson, Santa Clara, CA 


V. E. Hoggatt, Jr., proved that no three distinct Fibonacci numbers can be 
the lengths of the three sides of a triangle. (See page 85 of Fibonacci and 
Lucas Numbers, Houghton Mifflin Mathematics Enrichment Series, Houghton 
Mifflin, Boston, 1969.) Since a Pythagorean triple gives integral lengths for 
the sides of a right triangle, his result is more general. Hoggatt's elegant 
proof follows, where a, b, and e@ are the sides of the triangle: 


In any triangle, we must have a+b>e, b+e>a, andet+a> b. 
For any three consecutive Fibonacci numbers, F, + F,4, = Fy,42, and so 
there can be no triangle with sides having measures F,, F,4 1, Fyx42- In 
general, consider Fibonacci numbers, F,, F,, F,, where F, < K-j and 
fo49 Sh, Since Fo. +h, = fa.) end 1, S Pas We Rave #2 Shans 
and since F,,, < F;, we have F, + Fs < Fy. Therefore, there can be no 
triangle with sides having measure F,, F,, and F,. 


Also solved by Charles Ashbacher, Paul S. Bruckman, Piero Filipponi, C. 
Georghiou, Sahib Singh, Lawrence Somer, and the proposer. 


Each Term a Multiple of 3 
B-611 Proposed by Herta T. Freitag, Roanoke, VA 


Let 
n 


S(n) = Lupe? 
k=1 
For which positive integers n is S(m) an integral multiple of 3? 
Solution by Bob Prielipp, U. of Wisconsin-Oshkosh 


We shall show that S(n) is an integral multkple of 3 for each positive 
integer n. 
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The claimed result is an immediate consequence of the following lemma. 


Lemma: 3 divides Ly,42 for each nonnegative integer k. 


Proof: Because Ly = 3, the specified result holds when kK = 0. Let j be a non- 
negative integer. Then 


Luc g+1y)42 = Lugee = Zugeu + Luzes 

= QLusaa + Luge) + Guzee + 2Ly 343) = 2Ly-49 + 3L 
Hence, if 3 divides Ly;,2, then 3 divides Ly (3+ )42- 
follows by mathematical induction. 


4j7+3° 
The required result now 


Also solved by Paul S. Bruckman, Piero Filipponi, C. Georghiou, L. Kuipers, 


Chris Long, Br. J. M. Mahon, H.-J. Seiffert, Sahib Singh, Lawrence Somer, 
H. J. M. Wijers, Gregory Wulczyn, and the proposer. 


When the Sum Is a Multiple of 7 
B-612 Proposed by Herta T. Freitag, Roanoke, VA 


Let 


n 


T(n) “i Puxs2° 
=1 


For which positive integers m is T(m) an integral multiple of 7? 
Solution by Lawrence Somer, Washington, D.C. 
By inspection, we observe that the period of {F,,} modulo 7 is 16. 


F, = 121 (mod 7), Fe = 8 = 1 (mod 7), 


Fig = 55 = -l (mod 7), Fux = 377 = -1 (mod 7). 
It thus follows that 


Now, 


Fupgo = 1 (mod 7) if k = 0 or 1 (mod 4) 
and 


Fungo 2-1 (mod 7) if k = 2 or 3 (mod 4). 


Consequently, it follows that T(n) is an integral multiple of 7 for a positive 
integer m if and only if nm is an even integer. 


Also solved by Paul S. Bruckman, Piero Filipponi, C. Georghiou, Br. J. M. 


Mahon, Bob Prielipp, H.-J. Seiffert, Sahib Singh, H. J. M. Wijers, Gregory 
Wulczyn, and the proposer. 


Finding the Constants 
B-613 Proposed by Piero Filipponi, Fond. U. Bordoni, Rome, Italy 
Show that there exist integers a, b, and ec such that 
Fe, + F2_, = akere + b(-1)’F2 + c(-1)"FB. 
Solution by C. Georghiou, University of Patras, Greece 


We will show that a = 5 and b = c@ = 2. Indeed, from the identity 
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5F2 = Lo, - 2(-1)", 
we find 


5F2. + 5F2 =2L +L 


n+p n-p 2n+2p 
272 = rey ig Sa rie 
25h FS Lo Hop 2(-1) Lo, 2(-1) Lo, + 4(¢ 1) ‘ 
It follows, therefore, that 


Fe, , + F2_, - SF2F2 


n+p _ ss -1)"tP 
2n-2p = 4(-1) er Lo, hop 4( 1) 
and 


(2(-1)?Z,, + 2(-1)"L,, - 8(-1)"*?)/5 
2(-1)?F2 + 2(-1)"Fo. 


P 


r 


Also solved by Paul S. Bruckman, Herta T. Freitag, L. Kuipers, Br. J. M. 
Mahon, Bob Prielipp, H.-J. Seiffert, Sahib Singh, Lawrence Somer, Gregory 
Wulczyn, and the proposer. 


Quadruple Products Mod 8 
B-614 Proposed by Piero Filipponi, Fond. U. Bordoni, Rome, Italy 


F Show that 


n+1 


P. 


Let L(n) = L, oly -yln4 yh ii 


n-2 n-l'- n+l 


L(n) = F(n) (mod 8) 


and express [L(n) - F(n)]/8 as a polynomial in F,,. 


nik and F(n) = Taine at 


Solution by Sahib Singh, Clarion U. of Pennsylvania, Clarion, PA 


Using I, and I,, in Hoggatt's Fibonacci and Lucas Numbers, we get: 
L(n) = L#- 25 and F(n) = FH- 1. 
Replacing L2 by 5F2 + 4(-1)", we get 
L(n) - F(n) = 24F4* + 40(-1)"F2 - 8 = 0 (mod 8). 


Hence, 


LG) - Fin) _ am -1)"p2 — 
5 = 3F4 + 5(-1)"F2 - 1. 


Also solved by Paul S. Bruckman, Herta T. Freitag, C. Georghiou, L. 
Kuipers, Bob Prielipp, H.-J. Seiffert, Gregory Wulczyn, David Zeitlin, and 
the proposer. 


Identity for Iterated Lucas Numbers 
B-615 Proposed by Michael Eisenstein, San Antonio, TX 


Let C(n) L, and a, = C(C(n)). For n=0, 1, ..., prove that 


Bnt3 ~ Fn42%41 + a,- 
Solution by C. Georghiou, University of Patras, Greece 


= gh 4 gli), 


It is easy to see that a, Therefore, 


(q2(n+2) - gL(nt2)) (glint) ip gl(n+1)) 


B+ 2% 41 
= gh +3) 4 gL (n+3) Es (-1)2*F D (gh af gh(")) 


from which the assertion follows. 
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Also solved by Paul S. Bruckman, Piero Filipponi, Herta T. Freitag, L. Kui- 
pers, Br. J. M. Mahon, Bob Prielipp, H.-J. Seiffert, Sahib Singh, Lawrence 
Somer, David Zeitlin, and the proposer. 


KKK 
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ADVANCED PROBLEMS AND SOLUTIONS 


Edited by 
Raymond E. Whitney 


Please send all communications concerning ADVANCED PROBLEMS AND 
SOLUTIONS to RAYMOND E. WHITNEY, MATHEMATICS DEPARTMENT, LOCK 
HAVEN UNIVERSITY, LOCK HAVEN, PA 17745. This department especially 
welcomes problems believed to be new or extending old results. Proposers 
should submit solutions or other information that will assist the editor. To 
facilitate their consideration, all solutions should be submitted on separate 
signed sheets within two months after publication of the problems. 


PROBLEMS PROPOSED IN THIS ISSUE 
H-428 Proposed by Larry Taylor, Rego Park, NY 


Let j, m, and m be integers. Let qa and b be relatively prime even-odd in- 
tegers with b not divisible by 5. Let A, = aL, + bF,. Then A, = Any, - An-] 
with initial values 4; = b+a, A_-; =b-a. 

Prove that the following three numbers 


Ft hgage Ty eh PAS 


n+j m+ j? 2j°ntm 


are in arithmetic progression. 


H-429 Proposed by John Turner, Hamilton, New Zealand 


Fibonacci enthusiasts know what happens when they add two adjacent numbers 
of a sequence and put the result next in line. 

Have they considered what happens if they put the results in the middle? 

They will get the following increasing sequence of T-sets (multi-sets): 


T, = {1} 
given initial sets 
f, = {1, 2} 
T'3 - ti, 3; ys 
h.4, 3, $3285 
Ts = th. aa i. Ps ae 8, 2 Ts ae 
+~_* fie ty. Sy By Me hs Ty BS TB i 28 Fs BD 
etc. 


Prove that for 3 < 7 < n the multiplicity of 7 in multi-set 7, is 46(7), where 
¢ is Euler's function. 


SOLUTIONS 
What's the Point? 
H-406 Proposed by R. A. Melter, Long Island U., Southampton, NY 


and I. Tomescu, U. of Bucharest, Romania 
(Vol. 25, no. 1, February 1987) 
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Let A, denote the set of points on the real line with coordinates 1, 2, 


-» n. If F(n) denotes the number of pairwise noncongruent subsets of A,, 
then prove 


2n-2 4+ 2"/2 ~ 1] for even n, 
F(n) = 
27-2 + 32 2¢"-3)/2 ~ 1 for odd n. 


Solution by the proposers 
Let m =n, + +++ +n, be a decomposition of m into k nonnegative parts. It 
is well known that the number of such decompositions is equal to 


(4) 


The decompositions , +--+ +m, and nm, + +++ +n, will be said to be conju- 
gate. 
It follows that 


n-1 
F(n) = 1+ } a(m) 
m=1 


where a(m) is the number of pairwise nonconjugate decompositions of m. 


Denote by a(m, k) the number of pairwise nonconjugate decompositions of m 
into k parts: 


We shall consider four cases: 
(i) k = 22 and m = 2p. 


In this case the number of self-conjugate decompositions of m with k parts 
is equal to 


oe 
and hence, 


GoGo ae ee 


nous 06m, 10 = (E23) ~ 32) - @)] = 3G) #20) 
In order to calculate a(m), consider two subcases. 


I. Let m= 2p. It follows that 


om PEC) + Bs) [ED Zo wl 


220 Q21 Q21 


AZ CH) +E ()] = pee? + 2) = 2te-8 + ae 


m-2 (m-2)/2 
2 -2 . 


II. Let m = 2p + 1. One can write 


sEln™ 2) +H EG) * BG) - FEF) + 2B) 


220 220 220 


a (m) 


gm-2 4 9(m-3)/2, 
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aha oe, 
“32 1) +30 - 1) 


(41) k= 22, m= 2p+1. 


a(m, k) 


Here there are no self-conjugate decompositions; hence, 


a(m, k) = 3(7 _ 3): 


(iii) k = 22 +1, m= 2p. 


In order to count the number of self-conjugate decompositions, observe that 
the central position (m,) of m, + +++ + Moo4, must be an even integer. 
Thus, the number of self-conjugate decompositions is equal to 


Co ae ee 
It follows that, in this case 
Gee ae 
Hn tH 4) 
(iv) K=22+1, m= 2p+1. 


a(m, k) 


It can be seen that the central position of a self-conjugate decomposition 
must be an odd number. ; 


Finally, for odd n, 
F(n) = 1 + 27242714204 204 214204 22421 4 00. + 2%-3 + 20%-3)/2 
eS EE” icc ey PEPE Ol 2 + 06s 
$+ 2in-5)/2 4 g(n-5)/2) 4 2(n-3)/2 
m 9n-2 4 30 2("-3)/2 ~ 1, 


For even n, one obtains 
F(n) = 1+ 272 + 271 + 29 + 29 + 21 4+ 29 + 22 + 21+... 
gna On) 2: aot =3) en Sie 
= 2+ (29 +... + 2773) + (29 + 20 + 00. + 2Um-H/2 4 QKn-4)/2) 
ee ee Oe eB 
Also solved by Paul Bruckman. 


Nice End Product 


H-407 Proposed by Paul S. Bruckman, Lynwood, WA 
(Vol. 25, no. 1, February 1987) 





Find a closed form for the infinite product: 


= (5n + 2)(5n + 3) 
wl (5n + 1)(5n + 4)° (1) 
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Solution by Carl Libis, student, Tempe, AZ 


Use Theorem 5, p. 14, in Rainville's Special Functions, which says: 


8s 
> 5,, and if no a, or b, is a negative integer, 


s 

" 

If 2, a, 
k=1 k=1 


=o (n+a,)(n+a,) ... (N+a,) TCL +b,))T(1 +b,) ... TL +5,) ,, 
n=l (n+b,)(n +b)... (n+b,) Tl +a,)r(l+a,) ... Tl +a,) 











Thus, a 9 3 
= (5n + 2)(5n + 3) _ 5°95 il (n + s)(n + 5) 
n=0 (5n + 1)(5n + 4) z 2 eat @ ‘ ?) - +) 
373, s)r(5) 
573 °(3)r(5) 
573 | ()rG)(5)r65) 
573 (5)°(5)(s)45) 
QE) se) 


Also solved by D. Antzoulakos, O. Brugia & P. Filipponi, C. Georghiou, W. 
Janous, B. Prielipp, J. Shallit, and the proposer. 


Ghost from the Pasi 


H-125 Proposed by Stanley Rabinowitz, Far Rockaway, NY 
(Vol. 5, no. 5, December 1967) 





Define a sequence of positive integers to be left-normal if given any 
string of digits, there exists a member of the given beginning with this string 
of digits, and define the sequence to be right-normal if there exists a member 
of the sequence ending with the string of digits. 

Show that the sequences whose mth terms are given by the following are 
left-normal but not right-normal. 


P(n), where P(x) is a polynomial function with integral coefficients. 
P,,» where P, is the nt® prime. 

n! 

F,, where F, is the nt® Fibonacci number. 


andap 


Comment by Chris Long, student, Rutgers U., New Brunswisk, NJ 


Left-normality for all of the above was established by Raymond E. Whitney 
in this journal, vol. 11, no. 1, p. 77, and vol. 11, no. 2, pp. 186-1873; he 
also established that (b) and (c) are not right-normal. For (d), note that the 
Fibonacci sequence is defective mod 8, and hence is defective mod 1000; this 
shows that the Fibonacci sequence is not right-normal. However, the statement 
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that no polynomials with integer coefficients are right-normal is false, as the 
example P(n) = n demonstrates. Indeed, David Moews, a student at Harvard 
College, came up with the following characterization of right-normal polyno- 
mials with integer coefficients. 


Theorem (Moews): If Q is a polynomial with integer coefficients, then Q is 
right-normal iff 

(1) for all m there exists n with Q(n) =m (mod 10); 

(2) for all n, (Q'(n), 10) = 1, where Q'(x) is the formal derivative of Q(x). 
Proof (Moews): Note that for all m 21, Q@ can be viewed as a function from 
Z/10O"Zinto Z/10"Z. Q will be right-normal just when this function is 
surjective for all m; since Z/10"Z is finite, this will be the case just when 


this function is injective for all m. We induce on m to show that (1) and (2) 
imply this. 


If m= 1, this is clear; otherwise, let m> 1. Suppose we have x, y with 
Q(x) = Q(y) (mod 10"). Then Q(x) = Q(y) (mod 10"-!), so by the induction hy- 
pothesis, x = y (mod 10"-1). Let x = y + k10"-!. Then, since m > 1, 2(m-1) 2 
m, so Q(x) = Q(y) + k10"-1Q'(y) (mod 10”), and we must have k10"-1Q'(y) = 0 (mod 
10”), i.e., kQ'(y) =0 (mod 10), which gives k = 0 (mod 10) since Q'(y) is rela- 
tively prime to 10. Hence, x = y (mod 10”), which completes the induction. 


For the other implication, it is clear that Q@ cannot be right-normal if (1) 
fails. If (2) fails, let n have (Q@'@), 10) = a, a>1. Then, 1£ 5 = 10/4; 


Q(n + 10b) = Q(n) + 10bQ'(n) (mod 100), 
and a divides Q'(n) so 100 divides 10bQ'(n), which means that 
Q(n + 10D) = Q(n) (mod 100). 


This proves that @ is not injective as a function from Z/100Z to Z/100Z, so Q 
cannot be right normal. Q.E.D. 


Examples that David Moews came up with of Q(x)'s which satisfy conditions 
(1) and (2), and that are therefore right-normal, include 


Q(x) = ax + b for (a, 10) = 1 
and higher degree polynomials such as 


Q(x) = 2x5 + Sat + 5x2 + On. 


KAAKAK 
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